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Summary 

The computer program Coolant Passage Flow (CPF) was 
developed specifically for radial turbomachinery but can be 
used to analyze any turbomachinery coolant flow path geom- 
etry that consists of a single flow passage with a single inlet 
and exit. Flow can be bled off for tip-cap impingement cooling, 
and a flow bypass can be specified in which coolant flow is 
taken off at one point in the flow channel and reintroduced 
at a point farther downstream in the same channel . The user 
may either choose the coolant flow rate or let the program 
determine the flow rate from specified inlet and exit conditions. 
The main coolant flow must be subsonic everywhere except 
at the exit, where choked flow is allowed. Choked flow is also 
allowed at the bypass exit and at the tip-cap inlet and exit. 

CPF integrates the one-dimensional momentum and energy 
equations along a defined flow path and calculates the coolant s 
flow rate, temperature, pressure, and velocity and the heat 
transfer coefficients along the passage. The equations account 
for area change, mass addition or subtraction, pumping, 
friction, and heat transfer. 

The flow path geometry is described in Cartesian coordinates 
at node points that, in turn, divide the flow path into a desired 
number of intervals. Each interval may be further subdivided 
into any desired number of slices to increase calculation accur- 
acy. Standard correlations are used to predict the frictional 
pressure drop and the heat transfer coefficient for each interval, 
which is specified as being either a plain, a pin-finned, a trip- 
strip, or a finned passage. Alternatively, the friction factors 
and the heat transfer coefficients for each interval may be either 
specified directly or input as curves in tabular form. 

This report presents a detailed description of the program 
input and output, including error messages. Three example 
problems illustrate the use of the code, demonstrate the 
calculation accuracy as a function of the change in free flow 
area between nodes, and show a comparison with Fanno and 
Rayleigh line predictions. The report also presents all the 
analysis equations, as well as the derivations of the one- 
dimensional momentum and energy equations. 

Introduction 

In small turbomachinery the radial turbine has several 
inherent advantages over an axial turbine. Among these are 
higher overall work extraction, greater aerodynamic efficien- 


cy, and less sensitivity to tip clearance. These attributes have 
resulted in the widespread use of uncooled radial turbines 
in automotive turbochargers and auxiliary power units, but 
their use in primary powerplants (where higher operating 
temperatures and pressures are desired) has been held back 
by the difficulty of fabricating a cooled radial turbine. Both 
the U.S. Army and NASA, in a series of contractual efforts 
(refs. 1 to 4), have investigated techniques for fabricating 
cooled radial turbines. These efforts, although contributing 
significantly to solving the manufacturing problems, focused 
relatively little effort on the internal heat transfer of cooled 
radial turbines. The internal heat transfer calculations were 
performed, for the most part, by using company proprietary 
flow codes, along with approximations based on axial turbine 
experience. When NASA and the Army Propulsion Directorate 
decided to design an in-house cooled radial turbine for possible 
evaluation in the new Small Warm Turbine Facility, a 
computer code had to be developed that would accurately 
predict coolant flow and heat transfer inside cooled 
turbomachinery. The computer code Coolant Passage Flow 
(CPF) described in this report is the result of that effort. 

Although specifically developed for radial turbomachinery , 
CPF can be used to analyze any turbomachinery coolant flow 
geometry that consists of a single flow passage with a single 
inlet and exit. CPF integrates the one-dimensional momentum 
and energy equations along a defined flow path and calculates 
the coolant’s flow rate, temperature, pressure, and velocity 
and the heat transfer coefficients along the passage. The 
equations account for area change, mass addition or subtrac- 
tion, pumping (radius change), friction, and heat addition. 
Flow can be bled off for tip-cap impingement cooling, and 
a flow bypass can be specified in which coolant flow is taken 
off at one point in the flow channel and reintroduced at a point 
farther downstream in the same channel. The user may either 
choose the coolant flow rate or let the program calculate the 
flow rate from specified inlet and exit conditions. 

This report describes the program’s input, output, and error 
messages and illustrates the use of the code by three example 
cases. All the analysis equations are shown, along with their 
derivations and necessary assumptions where appropriate. 
An added feature of CPF that allows direct specification of 
heat generation in a given flow interval is described in appen- 
dix E, written by Teresa R. Kline. A heat generation term has 
been added to the energy equation. This feature allows the 
direct solution of the Rayleigh line problem (see example 
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problem 3). CPF can be obtained from the Computer Software 
Management Information Center (COSMIC), University of 
Georgia, Athens, Georgia, 30602. 


Flow Path Modeling 

A composite sketch of the types of flow and heat transfer 
configurations that can be analyzed with CPF is shown in 
figure 1. CPF determines the coolant’s flow rate, temperature, 
pressure, and velocity and the heat transfer coefficients along 
any arbitrary passage with a single inlet and exit. Note that 
internal flow division is not allowed, except for the shown tip- 
cap impingement flow and bypass flow. Bypass flow is taken 
off at one point in the flow channel and reintroduced at a point 
farther downstream in the same channel. The three flow 
channels shown in figure 1 would thus be analyzed separately. 
The calculation frame of reference is the rotating frame 
(moving with the turbomachinery). Thus, when the rotational 
speed is greater than zero, all coolant input and output total 
pressures and temperatures are relative total conditions. 

The flow path is defined with the Cartesian coordinates z, 
r, and y as shown in figure 1, where z is the axial direction 
and the axis of rotation, r is the radial direction, and y is 
perpendicular to both z and r. (All symbols are defined in 
appendix A.) The flow path geometry is described at nodes 
(at the passage centroid) that, in turn, determine a specified 
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number of intervals (input parameter NI). The nodes and 
intervals for the example 1 upper flow path are shown in 
figure 2. Along the flow path, the interval is that flow 
path length immediately preceding node n. Note that the inlet 
is specified separately and that the last (NK*) node is the exit. 
Also note that several nodes are coincident. These coincident 
nodes result from specifying zero-length “fake” intervals. 
These intervals, along with all input variables, are fully 
described in the section Program Input. 

Flow path geometry variables consisting of node location, 
free flow area, and hydraulic diameter are specified at each 
node. Standard correlations are used to predict the frictional 
pressure drop and the heat transfer coefficient for each interval, 
which is specified as being either a plain, a pin-finned, a trip- 
strip, or a finned passage. Alternatively, the friction factors 
and the heat transfer coefficients for each interval may be either 
specified directly or input as curves in tabular form. Each 
interval must be designated either a regular or a fake interval. 
Sudden step changes in flow area are modeled via fake 
intervals, in which the total pressure drop is calculated from 
established correlations. 

The change in free flow area between nodes defining a 
regular interval has a strong influence on the accuracy of the 
flow calculations. This means that, for geometries with 
significant area change, the nodes should be spaced close 
together but, for geometries with small area change, the nodes 
may be spaced farther apart. In order to avoid having to input 
an excessive number of nodes, an option can be used that 
divides each input interval (the interval that lies between input 
nodes) into a specified number of slices to get a more accurate 
flow solution. The program then switches from the input nodal 
frame of reference (shown in fig. 2) to the new, slice frame 
of reference, and new node and new interval numbers are 
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Figure 2. — Nodal representation of Figure 1 upper flow path. 
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created in accordance with the total number of slices. All 
calculations are then performed in the new frame of reference. 
The required node parameters needed at the new' node (slice) 
locations are obtained by linear interpolation of the appropriate 
original input nodal values. Input variables that were specified 
for the original input intervals are retained for the new intervals 
(which were created by slicing up the original intervals). For 
additional calculation accuracy, the program further divides 
all intervals (either the original input intervals or the calculated 
intervals generated by the slice option) into smaller sections 
called steps. The flow solution is marched along these steps 
from node to node and iteratively converged on the flow Mach 
number. The iterative flow solution procedure is described in 
detail in appendix B. 

Output may be specified in the original frame of reference 
or in the new, slice frame of reference, The total number of 
slices allowed has arbitrarily been set at 250. This can be 
increased to any desired number by changing the DIMENSION 
statements at the beginning of the code. Example 2 (sec 
the section Example Problems) illustrates the calculation 
accuracy achieved for various rates of free flow 7 area change 
between nodes. 

The user may either choose the coolant flow rate or let the 
program determine the flow rate from specitied inlet and exit 
conditions. The coolant flow must be subsonic at all internal 
flow path nodes. Choked (sonic) flow is allowed only at the 
exit node, the bypass exit, and the tip-cap inlet and exit. If 
the coolant flow is calculated from specified inlet and exit 
pressures and supersonic flow is calculated at internal flow 
path nodes, the weight flow' is reduced to achieve subsonic 
flow'. If a weight flow is specified and that flow is too large 
for the flow geometry, the program terminates with an 
appropriate error message. 

Program Input 

The program input parameters are listed in appendix C, 
along with the parameter description, the required units, the 
type of variable (real or integer), and the default values. Except 
for the title, all data input is in NAMELIST format. Although 
the appendix is essentially self-explanatory, the following 
matters and variables require further discussion. 

Input Curves 

Various parameters may be specified as curves in tabular 
form. Each table is described by five input parameters: an 
integer variable describing the number of points in the tabic, 
two subscripted parameters corresponding to the curve x and 
v values, and two parameters that specify the type of x and 
v axes (Cartesian or logarithmic (base 10)). Subroutine 
SPLFIT generates a spline curve fit from the tabular input. 
Each curve must be described by a minimum of 3 to a 
maximum of 25 points, with the abscissa fv axis) values input 
in ascending order. The slopes at the curve end points are 


evaluated from the first two and last tw r o data points. These 
points must thus be selected to give a good representation of 
the curve slope at each end. 

Flow Angles 

Figure 3 is a schematic of the angle definitions for general- 
ized bypass and impingement flows. At present the angle inputs 
ALIMP and ALBPI are dummy variables that have no effect 
on any calculations. They have been included for anticipated 
future uses as explained in the description of those variables. 

Nonsubscripted Variables 

Tide.— The title input may be up to 80 alphanumeric charac- 
ters in length. It is the overall title for each job submission 
and does not get input repeatedly if successive cases are run 
in a single data stack. The title is, however, printed out tor 
each case in a data stack. 

IOPPRP .— If the physical properties specification is set 
equal to 0, the user inputs the coolant's specific heat ratio, 
specific heat, viscosity, and thermal conductivity. These values 
are then assumed to remain constant for all flow calculations. 
If IOPPRP is set equal to 1 , the coolant's physical properties 
are evaluated from the subroutine AIRPRP. The properties 
vary and are calculated at local static temperatures. 

PTSPand TTSP — The inlet supply total pressure and total 
temperature must be specitied in the relative frame of reference 
(relative total conditions) if the rotational speed (input 
parameter RPM) is greater than 0. For no rotation these 
variables are specified in the absolute frame of reference. 
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Figure 3.— Angle definiiions lor generalized How path. 
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KTIN .— The inlet total pressure loss coefficient accounts 
for various inlet configurations and can assume a value 
between 0.0 and 1 .0. If it is equal to 0.0, the inlet total 
pressure is set equal to the specified supply total pressure 
(PTSP); if it is equal to 1 .0, the inlet total pressure loses one 
full velocity head. 

Nodal or Interval Variables 

Geometry parameters or geometry and flow and heat transfer 
parameters must be supplied at the input nodal points and for 
the intervals defined by these nodal points. For the nodal or 
interval subscripted parameters of appendix C (identified as 
variable type R(NI) or I(NI)), only the variables Z, R, Y, A, 
and DH are nodal values, in that they define values at each 
input node location. All other variables are interval values and 
represent average values over the interval. The nodal values 
in appendix C are identified by the superscript nd , and the 
interval values by the superscript /. 

For each nodal or interval input parameter, an array of NI 
numbers must be specified, even though actual values may 
be needed for only a few nodes or intervals. The values that 
are not needed are thus dummy or fill-in values, but they must 
still be specified. For example, if a flow path consists of 20 
intervals and the user wishes to let the program calculate the 
friction factors in the first 10 intervals but wants to specify 
the friction factors in the last 10 intervals (via the parameter F), 
the first 10 numbers of the 20-number F array are dummy 
values, and only the last 10 numbers are used by the code. 

NSLI(n ).— For large area variations between nodes the 
calculation procedure introduces small artificial increases or 
decreases in calculated total pressure and total temperature (see 
example 2 in the section Example Problems). This parameter 
specifies the number of slices into which each interval is to 
be divided in order to increase the accuracy of the calculation 
procedure. For certain intervals, only one slice (no dividing) 
is allowed, and no matter what values of NSLI(/?) are input, 
the program defaults NSLI(n) to 1 . These intervals are (1) fake 
intervals; (2) tip-cap flow intervals; (3) bypass intervals (flow 
leaving and entering); and (4) intervals in which a total 
pressure loss coefficient ( K factor; K T ) is specified by the 
input variable KTSGMT. 

IOPTFI(n ).~ This parameter designates fake intervals, 
which are used to model abrupt area changes (such as occur 
when there is a change in the interval heat transfer configur- 
ation). By using fake intervals the correct total pressure drop 
can be calculated from established correlations for sudden 
increases or decreases in flow- area. Fake intervals may be used 
anywhere in the flow path and should be of zero length. That 
is, the coordinates of the nodes bracketing the fake interval 
should be identical. Example 1 (in the section Example 
Problems) shows how' fake intervals are used between intervals 
with different heat transfer geometries to model the sudden 
area changes. 

IOPTHT(n) and IOPTFF(n ). — These variables specify the 
types of heat transfer and friction factor calculations. 


respectively, for each interval. IOPTHT(m) and IOPTFF(n) 
are the main controlling parameters in that they determine the 
specific information used in each interval for the requested 
heat transfer and friction calculations. Only certain 
combinations of IOPTHT(n) and IOPTFF(*) are allowed for 
each interval. These combinations are shown in table I. 
The heat transfer coefficient may be user specified 
(IOPTHT(/?) — 0), obtained from internal correlations 
(IOPTHT (n) = 1, 2, 3, or 4 for trip strips, pin fins, finned 
passage, and plain passage, respectively), or may be deter- 
mined from an input Colburn 7- factor curve (IOPTHT(n) = 5). 
The friction factor may also be user specified (IOPTFF(n) = 0), 
obtained from trip-strip geometry input (IOPTFF(n) = 1), 
obtained from an appropriate input curve (IOPTFF(n) = 2), 
obtained from a roughened-plain-passage correlation 
(IOPTFF(n) = 3), or obtained from a specified K factor by 
means of the input parameter KTSGMT (IOPTFF(n) = 4). 
Only the combinations of IOPTHT(n) and IOPTFF(n) shown 
in table I are allowed. If unallowed combinations are chosen, 
the program aborts and prints out an appropriate error 
message. For fake intervals, dummy values must be input for 
both IOPTHT(n) and IOPTFF(n). 

Note that CPF does not check that the interval inputs are 
consistent with the specified values of IOPTHT(n) and 
IOPTFF(n). The user must thus make sure that the correct 
type of information is supplied for each input interval. 

FCOR( n), HCCOR( n), and HTACOR{n). —The user may 
modify values of friction factor, heat transfer coefficient, and 


TABLE I — ALLOWED COMBINATIONS OF 
IOPTHT AND IOPTFF 


Heat transfer 
coefficient option, 
IOPTHT 

Friction factor 
option. 
IOPTFF 

0— User specified 

0— User specified 

1— Trip strips 

1— Trip strips 

2— Pin fins 

0 — User specified 

2— Input curve (NPFF2, etc.) 

3 — Finned passage 

0 — User specified 

2— Input curve (NPFF3, etc.) 

3— Roughened-plain-passage 

correlation (ref. 7) 

4— Plain passage 

(internal correlation) 

f 

0 — User specified 

2— Input curve (NPFF4, etc.) 

3— Roughened-plain-passage 

correlation (ref. 7) 

4— Specified K factor 

5— Input Colburn /-factor 
curve (NPCJF, etc.) 

0— User specified 

2— Input curve (NPFF4, etc.) 

3— Roughened-plain-passage 

correlation (ref. 7) 

4— Specified K factor 
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heat transfer area by the correction factors FCOR(«), 
HCCOR(n), and HTACOR(n), respectively. Consider the flow 
path shown in figure 4, in which nodes n — 1 and n define 
a 180° bend. The program calculates the straight-line distance 
between the nodes. Note that the true distance between the 
nodes (along the flow direction) is significantly greater than 
the calculated distance. If the friction pressure drop is specified 
by the Fanning friction factor (IOPTFF (n) = 0, 1, 2, or 3), 
the friction factor must be modified by 

True distance 

FCOR(n) ; 

Calculated straight-line distance 

If the friction pressure drop is specified by a total pressure 
loss coefficient (IOPTFF(n) = 4), no correction is needed, 
since the loss coefficient is specified for a given geometry 
(e.g., a 180° bend) and thus is independent of actual flow path 
length. For the geometry shown in figure 4, the heat transfer 
area must similarly be modified by 

True distance 

HTACOR(fl) = : 

Calculated straight-line distance 

The heat transfer coefficient correction factor HCCOR allows 
the user to account for inlet or turning effects. The use of the 
described modifiers allows flexibility in applying the program. 
Note, however, that their use cannot be hidden, since they 
are always shown in the program output. 

IOPTCR(n).— This input parameter controls the option for 
calculating a heat transfer coefficient correction due to rotation. 
If IOPTCR(n) is set equal to 0 (default value), no correction 
is made, and the heat transfer coefficient rotation modifier 
HCCORC(n) (not an input; see the section Program Output) 
is defaulted to 1.0. If IOPTCR(n) is set equal to 1, the 
calculated or specified heat transfer coefficient is modified as 
described in appendix B. Note that the correction is valid only 
for purely radial flow in plain passages and for fluid 
temperatures lower than the wall temperature (the fluid cooling 
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n - 1 AND n 

Figure 4.— Comparison of true and calculated distances in 180° bend. 


the wall). The calculated correction factor HCCORC(n) may 
be greater or less than 1 .0 and is shown in the program output. 

EDD(n) and PDE( n). —These parameters describe the 
geometry for trip strips. They are shown schematically in 
figure 5. 

NFNP(n), THFN(n), SPFN( n), and HHTFNfa ).- These 
parameters describe the geometry for finnea passages. Tney 
are shown schematically in figure 6. 


INPUT PARAMETER 


DEFINITION 


EDD 

PDE 


e/D 

p/e 



INPUT PARAMETER 

DEFINITION 

THEN 

T f 

SPFN 

if + rn 

HHTFN 

Lf 

NFNP 

NUMBER OF FINNED 


PASSAGES 

i i 


/-plane through node n 

/ 


/ ^PLANE THROUGH NODE n + 1 
x / HEAT 

/ / / f FL LJX 
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A(n) and DH( n).— The flow area and the hydraulic diam- 
eter based on minimum free flow area must be input for each 
node. These nodal values are used in evaluating the local Mach 
and Reynolds numbers. 

AV( n) andDHV(n) .—The flow area and hydraulic diameter 
based on open volume are needed only for pin-fin intervals. 
They are required in addition to the input variables A (n) and 
DH(/?). The interval values AV(«) and DHV(n) are used in 
the pin-fin heat transfer correlation. A V(n) is equal to VI t 
where V is the open volume in the interval and f is the interval 
length; DH V(«) is equal to 4 K/S, where S is the total heat 
transfer surface area in the interval (endwall plus pin fins). 

Parameters Specifying Location of Tip-Cap 
and Bypass Flow 

The variables NTCF, NBPI, and NBPO define the interval 
in which a tip-cap flow, the start of a bypass channel, or the 
end of a bypass channel, respectively, is specified. This is 
illustrated in the example flow path of figure 2, where the tip- 
cap flow leaves from interval 10 and the bypass channel starts 
in interval 2 and ends in interval 18. The first interval should 
not be a tip-cap interval or contain the start of a bypass channel 
(NTCF ^ 1, NBPI ^ 1). Any one interval can only be 
assigned one of these three options. This means that a tip cap 
and the start of a bypass channel cannot be specified for the 
same interval (NTCF ^ NBPI), nor can a tip cap be in the 
interval ending a bypass channel (NTCF ^ NBPO). It is also 
not recommended that the last interval contain the end of the 
bypass channel (NBPO s* NI). Note that tip-cap or bypass 
flows cannot be specified for fake intervals. 

Further comments are in order regarding the use of a bypass 
channel. Subroutine BYPSFL, which calculates the bypass 
flow, is a simplified version of the main routine used to 
calculate the coolant flow in CPF. The two main simpli- 
fications are that heat transfer is not considered for the bypass 
flow and that all calculations are based on constant physical 
properties. Both the inlet and exit conditions for the bypass 
channel are taken to be the average of the values at the 
respective interval end points. This implies that, for maximum 
accuracy, these intervals should be kept as short as possible. 
Further accuracy w ill be gained if these intervals are of constant 
cross-sectional area. Also, for best results, the bypass flow 
should not be a significant fraction of the main coolant flow. 

Tip-Cap Parameters for Interval NTCF 

The parameters that describe the geometry for tip-cap flow 
are shown schematically in figure 7. The impingement angle 
ALIMP shown in figure 3 is a dummy variable that has no 
effect on the calculations. It is included for possible future use 
in a more general impingement flow calculation procedure. 

Bypass Channel Parameters 

PBPRI.— This variable specifies the pressure recovery at 
the bypass inlet and is determined by the geometry of that 
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Figure 7.— Tip-cap geometry parameters and flow stations. 


inlet. PBPRI must be set between 0.0 and 1 .0. Setting PBPRI 
equal to 0.0 causes the bypass inlet total pressure to be set 
equal to the bypass inlet interval average static pressure. 
Such a condition exists if the bypass flow is drawn off perpen- 
dicular to the main coolant flow. Setting PBPRI equal to 1 .0 
causes the bypass inlet total pressure to be set equal to the 
bypass inlet interval average total pressure. Such a condition 
exists if the bypass flow is drawn off tangential to the main 
coolant flow. 

ALBPL— This angle describes the component of velocity 
of the bypass exiting flow in the main-stream coolant flow 
direction for interval NBPI. ALIBPI is a dummy variable that 
has no effect on any calculation. It is anticipated that in the 
future this variable, along with additional geometry param- 
eters, will be used in a generalized calculation procedure for 
the pressure recovery at the bypass inlet. 

ALBPO.— This angle describes the component of velocity 
of the bypass entering flow in the main-stream coolant flow 
direction for interval NBPO. ALBPO is equal to 0° for bypass 
flow entering tangential to the main-stream flow and equal to 
90° for bypass flow entering perpendicular to the main-stream 
flow. Figure 3 shows all flow r angle definitions. 

Miscellaneous Comments 

During the initial flow path modeling the specified free flow 
areas, as well as the inlet and exit pressures, may not be 
compatible with physically possible flow . For complex flow 
geometries it is suggested that the initial calculations be 
performed with a specified flow rate (IOPTW = 0) to make 
the calculations independent of exit backpressure. Starting with 
a low flow rate (to ensure subsonic flow at all internal nodes), 
the high Mach number locations can be identified and the flow 
areas increased, if necessary, as the flow* rate is increased. 
Also, the calculated trailing-edge pressure can be used to size 
the exit area for compatibility with the backpressure. 

Some care must be taken with the last interval for choked 
exit flow. If the flow area is constant in the last interval 
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DLBP = 0°6 5*! IB DHBP=0 N 0757 1 FFBP = 0 . 10 , AINBP = .001, PBPKI-0.50, 
R1NBP=3 80, A0UTBP= .001, ROUTBP=3.75, ALBPO-36.U. 

10PHTT = i, NHI=2, DH^.030, NHO = 3, DHO=.020, PBP = 200.0, 


50, 

30, 


43. 

.03, 


2.35, 2.30, 126*0.0, 


. 0100 , 
.0300, 
AV=17*0 . 0, 
DH= . 2310 , 

. 0696, 

. 0888 , 


KMET = 23*Q . 0, 
0300, .0175, 

.0225, .0275, 

126*0 .0, 


12.0, 126*0.0, 
.0126, .0126, 
0500, .0300, 


0085, 
0300 , 


. 0085, 
.0300, 


. 1020 , 
. 0824, 
. 0980, 


.0820, 
.0857, 
. 0833, 


,0775, .0775. 
.0880, .0960, 
126*0,0, 


. 0696, 
.0823 . 


R 1H = 6 . 60 , R1PL=6.75, SPIH=0.2O 
DIC= .075, A1MP= .0589, TWIMP=1500 


AL I IIP = 9 0 . 
8 END 


Figure 8.— Example 1 input data. 


( A (NI - 1) = A(NI)) and the user specifies zero friction and 
zero heat transfer (no heat addition), the solution may not 
converge. Without friction or heat addition in the last interval 
(both of which drive the flow toward choking) the flow will 
also want to choke at node NI — 1, but this cannot be 
accommodated. The problem can be circumvented by 
specifying the exit free flow area A(NI) to be slightly smaller 
than A(NI - 1). 

A complete data input set (for example 1 in the section 
Example Problems) is shown in figure 8. This input describes 
the upper flow path of figure 1 and illustrates the use of 
both tip-cap and bypass flow. Details concerning the 
preparation of this data input set are given in the section 
Example Problems. 


Program Output 

The program output is described in figure 9, which shows 
the output for example 1 . The title is shown first, followed 
by statements concerning the chosen options for physical 
properties, coolant flow calculation, rotational speed, whether 
tip-cap or bypass flow was specified, and if so, in which input 
intervals. The next statement tells the user whether the output 
is for the original input frame of reference or for the revised 
frame of reference (generated from the slices). This statement 
is followed by a table that shows the specified nodes, the 
number of slices in each interval, and the slice number (new 
node number) corresponding to each original input node. This 
table is important, since all error messages refer to the slice 
frame of reference. Any problem leading to warning or abort 
messages in the new frame of reference can be traced back 
to the input frame of reference with this table. 

The next output is a listing of all input data curves, followed 


by any warning or abort messages and a statement quanti- 
fying the overall number of iterations needed to achieve flow 
convergence. 

A tabulation of mostly geometric parameters is next printed 
out for the following node or interval items: 

I node or interval number 

IOPTFI input option for fake interval specification 
IOPTHT input option for heat transfer specification 
IOPTFF input option for friction factor specification 
IOPTCR input option for heat transfer rotation correction 
specification 

Z specified or calculated (for NSLI > 1) node z value 

R specified or calculated (for NSLI > 1) node r value 

Y specified or calculated (for NSLI > 1 ) node y value 

A specified or calculated (for NSLI > 1) node free 

flow area 

DH specified or calculated (for NSLI > 1) node 

hydraulic diameter 

TW specified interval wall temperature 

DX calculated straight-line interval distance between 

nodes (for all fake intervals DX is set to 0.0) 

XTOT flow path length up to that node (summation of DX) 

KT specified total pressure loss coefficient KTSGMT 

for regular intervals if IOPTFF(I) = 4; or 
calculated total pressure loss coefficient for fake 
intervals if IOPTFI(I) = 1; if a friction factor is 
specified or calculated, KT is shown as zero. 

NSLI number of specified slices in an interval it output 
is in original input frame of reference; set equal 
to 1 if output is in new, slice frame of reference. 
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ORiGiS^AL PAGE IS 
OF POOR QUALITY 


EXAMPLE PROBLEM 1 (FIGURE 1 UPPER FLOW PATH) 

ENGLISH UNITS 

PHYSICAL PROPERTIES (G, CP, MU AND K) OBTAINED FROM SUBROUTINE AIRPRP 

COOLANT FLOW RATE DETERMINED FROM SPECIFIED INLET AND EXIT PRESSURES 

ROTATIONAL SPEED IS 0.570E+05 RPM 

BYPASS CHANNEL SPECIFIED BETWEEN INTERVALS 2 AND 18 

TIP CAP FLOW SPECIFIED IN INTERVAL 10 

OUTPUT IN ORIGINAL NODE/ INTERVAL FRAME OF REFERENCE 


INPUT 

SLICES 

SLIC 

NODE 

IN 

NO. 

NO . 

INTERVAL 


1 

2 

2 

2 

1 

3 

3 

2 

5 

A 

2 

7 

5 

1 

8 

6 

2 

10 

7 

2 

12 

8 

2 

IA 

9 

1 

15 

10 

1 

16 

11 

1 

17 

12 

2 

19 

13 

2 

21 

IA 

2 

23 

15 

2 

25 

16 

2 

27 

17 

1 

28 

18 

1 

29 

19 

2 

31 

20 

2 

33 

21 

2 

35 

22 

2 

37 

23 

1 

38 

2A 

2 

AO 


CARTESIAN X-AXIS 
CARTESIAN Y-AXIS 

THE INPUT FANNING FRICTION FACTOR CURVE 


REA 
0 , 1000E105 
0 . 5000E+03 
0. 1D00E+0A 
0.2000EH3A 
0 . 2500E+GA 
0 . 3000E+OA 
0 . 3500E+OA 
0 .AOOOE+OA 
0 .5000E+OA 
0 . 1000E+05 
0. 1000E+06 
0. 100DE+07 


FFF A 
0. 170GE+00 
0 . 3AOOE-OI 
0 . 170DE-01 
0 .8500E-D2 
0 .70DOE-02 
0 . 77 00 E - 02 
0 . 1 OOOE-Ol 
0 . I 070E-01 
0 . 1000E-01 
0 .8000E-02 
0. A300E-O2 
0. 2900E-02 


FOR PLAIN PASSAGES IS 


LOGARITHMIC X 
LOGARITHMIC Y 
THE INPUT COLBURN 
REJCF 

0 . 1000E+03 
0 . 1 DOGE+OA 
0. 2D00E+0A 
0. 3000E+0A 
0 .6000E+OA 
0 . IOOOE+05 
0. I0D0E+06 
0 . 1 Q00E+07 


-AXIS 

-AXIS 

J-FACTOR CURVE 
CJF 

0. AAOOE-OI 
0 . A500E-02 
0 . 2300E-02 
0 . 2600E-02 
0. 3200E-02 
0. 2900E-02 
0. I 900E-02 
0 . 130DE-02 


IS 


WARNING - IMPINGEMENT FLOW HAS NOT CONVERGED IN SUBROUTINE TIPIMP IN OVERAU FLOW ITERATION < 
WARNING - IMPINGEMENT FlOH HAS NOT CONVERGED IN SUBROUTINE TIPIMP IN OVERAU FLOW ITERATION 9 


CONVERGENCE ACHIEVED IN 27 OVERALL ITERATIONS 


I 


2 

3 

A 


6 

7 

8 
9 

10 

11 

12 

13 

1A 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2A 


OPTIONS ( I OPT ) l 

FI HT FF CR (IN) 


R 

(IN) 


r 

(IN) 


0 A 
0 A 
0 A 

0 A 

1 7 
0 I 
0 I 
0 1 
0 A 
0 5 

0 A 
0 1 
0 1 
0 1 
0 1 
0 1 

1 7 
0 2 
0 2 
0 2 
0 2 
0 2 
1 7 
0 3 


2 

2 

2 

2 

7 

1 

1 

1 

A 

A 

A 

1 

1 

1 

l 

1 

7 

0 

0 

0 

0 

0 

7 

3 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

D 

0 

0 

0 

0 

a 

0 

0 


0 . OOOOE+OO 
0 . AOOOE + OO 
0. 5000E+00 
0 . AAOOE+GO 
0. 350DE+00 
0 . 3500E+00 
0 .2300E+00 
0 .200QE+00 
0 .2000E+00 
0 .2200E+00 
0. 980QE+00 
0. 50QOE+G0 
0. 5DQQE+0D 
0. 5200E+00 
0 . 7200E+00 
0 . 9500E+00 
0 . 1A00E + 01 
0 . IA00E+01 
0 . 1630E+01 
0. 1920E+Q1 
0 . 2 550 E + 0 1 
0 . 3AQ0E+O1 
0 . A300E + 01 
0 . A300E + 01 
Q . 4620E+ 0 1 


0. 3850E+D1 
0 . 3850E + 0I 
0 . AOOOE+Oi 
0. A550E+01 
0. A650E + 01 
0. A650E+01 
0. 5350E+0I 
0.6020E+01 
0 .6A30E+01 
0 .6A30E + 0I 
0 -6A30E + 01 
0 .6A30E+01 
0 .6100E + 01 
0 . 5550E + 0] 
0. A950E+01 
0 . A550E+01 
0 . AGOOE+Ol 
0. AOOOE+OI 
0 . 3850E + 01 
0 -3600E + 01 
0.3250E+01 
0.30606+01 
0. 3030E+01 
0. 3030E+01 
0 . 3D3GE+Q1 


0. OOOOE+OO 
0. OOOOE+OO 

o . ooaoE+oo 

0 . OQQOE+OQ 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0. OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE+OO 
0 .OOOOE+OO 
0 . OOOOE+OO 
0 . 7 OOOE-0 1 
0 . 7000E-01 
0 . 1000E+00 
0. 1800E+00 
0 . A700E+00 
0. 1150E+01 
0 . 2350E+01 
0 . 235DE+0 1 


A 

( IN**2) 

0 . 625QE-01 
0 . 7 5O0E-01 
0.63QQE-01 
0.5AOOE-OI 
0. A500E-01 
0.3000E-01 
0 . 1750E-01 
0 . 1260E-01 
0.1260E-01 
0 . 8500E-Q2 
0 . 850QE-02 
0 . 1000E-01 
0. 1000E-D1 
0. UOOE-Ol 
0 .1750E-01 
0.2250E-01 
0 . 2750E-01 
0. 3000E-01 
0. 3000E-01 
0.3000E-01 
0 . 3000E-01 
0.3000E-01 
0.2600E-01 
0. 2000E-01 
0 .20006-01 


0.2310E+00 
0.2310E+00 
0 .1590E+00 
0 . 1560E+00 
0. 1530E+00 
0. 1020E+00 
0 .8200E-01 
0.7750E-01 
0 . 7750E-01 
0.696QE-Q1 
0.6960E-01 
0.7270E-01 
0.7270E-01 
0.7A60E-0I 
0.82A0E-D1 
0 .8570E-01 
0.S8D0E-01 
0.9600E-01 
0.8880E-01 
0 . 8880E-01 
0 . 8280E-01 
O.8U0E-01 
0 .6380E-OI 
0.9800 E-0 1 
0 . 8330E-0 1 


TW 
< F ) 


0. 1AOOE+OA 
0. 1AOOE+OA 
0. 1AOOE+OA 
0 . 1 AOOE+OA 
0 . 1 AOOE+O A 
0 . 1 AOOE + OA 
0 . 1 AOOE+OA 
0. 1AOOE+OA 
0 . 1 AOOE+OA 
0 . 1 AOOE + OA 
0 . 1 AOOE+OA 
0 . 1 AOOE + OA 
0 . 1 AOOE + O A 
0. 1A00E + 0A 
0. 1AOOE+OA 
0. 1AOOE+OA 
0. 1AOOE + OA 
0 . 1 AOOE + OA 
0 . 1 AOOE+OA 
0 .1AOOE+OA 
0 . 1AOOE + OA 
0 . 1A00E+0A 
0. 1AOOE+OA 
0 . 1AOOE+OA 


DX 

(IN) 


0 . A OOOE+OO 
0 . 1803E+00 
0 . 3551 E+00 
0 . 3132E+00 
0 . OOOOE+OO 
0 . 7 102E+00 
D.67Q7E+00 
0. A100E+00 
0.2000E-01 
0 . 2600E+00 
0.200DE-01 
0 . 3300E+ 00 
0 . 550AE+00 
0 . 6 3256+00 
0. A61AE+00 
0. 71AIE+00 
0. OOOOE+OO 
0 . 276ZE+00 
0 . 3912E+00 
0 . 7769E+00 
0.1105E+01 
0. 15D0E+01 
0 . OOOOE+OO 
0 . 5522E+00 


XTOT 

(IN) 


0 . AOOOE+OO 
0 .5803E+00 
0 . 935AE+00 
0 . 12A9E + 01 
0. 12A9E + 01 
0. 1959E+01 
0 .2629E+01 
0 . 3039E+ 0 1 
0 .3059E+01 
0.3319E+01 
0. 3339E+01 
0 . 3669E+0 1 
0. A220E+01 
0. A852E+01 
0 . 531AE+01 
0 . 6 028 E+ 0 1 
0 .6Q28E+01 
0 .630AE+01 
0 . 6695E+01 
0 . 7A72E+01 
0.8577E+01 
0. 1008E+02 
0. 1008E+02 
0. 1063E+02 


KT NSLI 


0 . 50Q0E+ 00 
0. OOOOE+OO 2 
0, OOOOE+OO 1 
0. OOOOE+OO 2 
0. OOOOE+OO 2 
0.1037 E+00 1 

0. OOOOE+OO 2 
0. OOOOE + OO 2 
0. OOOOE+OO 2 
0.16 50E + 01 1 

0. I650E+01 l 
0.1650E+01 1 

0. OOOOE+OO 2 
0. OOOOE + OO 2 
0. OOOOE+OO 2 
0. OOOOE+OO 2 
0. OOOOE+OO 2 
Q.69A4E-02 I 
0. OOOOE+OO 1 
0. DOQDE+OO 2 
0. OOOOE+OO 2 
0. OOOOE+OO 2 
0. OOOOE + OO 2 
Q.923QE-01 1 

0. OOOOE+OO 2 


Figure 9.— Example 1 program output. 
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TIP CAP FLOW CALCULATIONS IN INTERVAL OR SLICE 10 ( 3 ITERATIONS) 

FLOW RATE= 0.1860ZE-02 (LBM/S) 


2 

0.30000E-01 (IN) 
0.90000E+02 (DEG) 
0 . 31 1G9E+00 
0 . 8 A 38 8 E + 0 0 


NO. OF INLET HOLES = 

INLET HOLE DIAMETER = 

ALIMP 

INLET MACH NUMBER - 
INLET DISCH. COEFF. = 

CHAMB. TOT. PRESS. (3) = 
CHAMB. MACH NUMBR . <3)= 
CHAMB. TOT. PRESS. (4) = 
CHAMB. TOT. TEMP. (4) = 
IMPINGEMENT H (4) = 

NO. OF OUTLET HOLES = 3 

OUTLET HOLE DIAMETER 3 0 

OUTLET MACH NUMBER = 0 

OUTLET DISCH. COEFF. = 0 

OUTLET BACK PRESS. - 0 


0 , 22927E+03 (PSIA) 

0 . 1920ZE-0 1 
0.23737E+03 (PSIA) 
0.13S85E+04 (F) 

0 14289E+04 ( BTU/FT»*2/HR/R ) 


. 20000E-01 (IN) 
.51 1 12E+00 
. 891 1 1E+00 
.2G0QOE+03 (PSIA) 


I 

T 

TT 


(F) 

(F) 


0 .7497E+03 

0 .7500E+03 

1 

0 .7700E+03 

0.7701E+03 

2 

0 .7765E+03 

0 .7766E+03 

3 

0 .8059E+03 

0 .8Q61E + 03 

4 

0 .830ZE+03 

0 .8305E+03 

5 

0.8298E+03 

0 .8305E + 03 

6 

0 . 1078E+04 

0 . 1 080E+04 

7 

0. 1212E+04 

0. 1216E+04 

8 

0 . 1264E+04 

0 . 1 268E+04 

9 

0. 1259E+04 

0 . 1268E+04 

ID 

0 . 1 26 3E + 04 

0. 1269E+04 

n 

0. 1265E+04 

0 . 1270E+04 

12 

0. 1296E+04 

0 . 13Q2E+04 

13 

0.1331E+04 

0 . 1339E+04 

1 4 

0 . 1360E+04 

0 . 1365E+04 

15 

0 . 1372E+04 

0 . 1377E+04 

16 

0 . 1353E+04 

0 . 1 388E+04 

17 

0. 1384E+D4 

0.1388E+04 

18 

0. 1257E+04 

0 . 1262E+04 

19 

0. 1274E+04 

0 . 1280E+04 

20 

Q.1302E+04 

O. 1 309E+04 

21 

0 . 1329E + 04 

0.1 339E+04 

22 

0 . 1342E+04 

0. 1365E+04 

23 

0.1325E+04 

0 . 1365E+04 

24 

0 . 1331E+04 

0 . 1372E+04 


P 

(PSIA) 

0.1748E+03 
0 . 1748E+05 
0 . 1800E + 03 
0 . 1932E+03 
0 . 2Q60E+03 
0. 2Q57E+03 
0.2310E+03 
0.2471E+03 
0 . 2549E+03 
0 . 24 59 E + 0 3 
0 .2432E + 03 
0 .2394E+03 
0 . Z084E+03 
0 . 1610E+03 
0 . 1256E + 03 
0.1101E+03 
0 . 9268E+02 
0 . 9282E + 02 
0 .9076E + 02 
D . 8595E+02 
0 .7827E + 02 
0 . 697 1 E + 02 
0 . 5320E+02 
0 . 509BE+02 
0 .5000E+02 


PT 

(PSIA) 

0 . 1749E+03 
0 . 1749E+03 
0.1801E+03 
0 . 1933E+03 
0.2061E+03 
0 . 206 1E + 03 
0 . 2322E+03 
0 . 2494E+03 
0 .2572E+03 
0 . 25L2E+Q3 
0 . 2464E+03 
0 . 2418E + 05 
0 . 21 12E+03 
0 . 164LE+03 
0 . L274E+03 
0. 1112E + 03 
0 . 9355E + 02 
0 .9554E+02 
0 .9185E+02 
0 .8711E + 02 
0 .7957E+02 
0.7119E+02 
0 . 5583E+02 
0 . 5564E+02 
0 .5478E+02 


M 

(LBM/5) 

0 .9745E-02 
0 .97A5E-02 
0 .5140E-02 
0 .8140E-Q2 
0.8140E-02 
0.8I90E-02 
0. 8140E-02 
0.8140E-02 
0 .8140E-02 
0 .8140E-02 
0 .6280E-02 
0 .628DE-D2 
0 .6280E-02 
0 .6280E-02 
0.6280E-02 
0.6280E-02 
0.6280E-02 
0 .6280E-02 
0 .7885E-02 
0 .7885E-02 
0 .7885E-02 
0 .7885E-02 
0 . 7885E-02 
0.7885E-02 
0 . 7885E-Q2 


V 

(FT/S) 

0.5757E+02 
0 .4877E+02 
0 .4735E+02 
0 .5269E+02 
0 .6045E+02 
0 .9077 E+02 
0 . 1652E+03 
0.2333E+Q3 
0 . 233LE+03 
0 . 3572E+03 
0 .2793E + 03 
0 .2414E + 03 
0 .2823E+0J 
0 .3387E+D3 
0 .2768E+03 
0 . 2478E+03 
0 .24236 + 03 
0 . 2219E+Q3 
0 . 2653E+03 
0 . 2830E+03 
0.3156E+Q3 
0 .3599E+03 
0 .5482E+03 
0 .7364E+03 
0 .7533E+03 


M 


0.34L6E-01 
0.2871E-01 
0.2780E-01 
0.3060E-01 
0.3979E-01 
0.52Z4E-01 
0 .8753E-0I 
0 . 1188E+Q0 
0. 1171E+0D 
0. 1796E+00 
0.1403E+00 
0.1212E+00 
0 . 1 405E+00 
0.1671E+00 
0. 1355E + 00 
0 . 1209E + 00 
0.1179E+00 
0 . 1079E+00 
0. 1335E+00 
0 . 1417E+00 
0.156 9 E + 00 
0.1776E+00 
0 . 2695E+00 
0 .363BE+00 
0.3716E+00 


GAMMA CP VISC. TH. LNDI. 

BTUZ(LBM-R) IBM/(FT-S) BTU/CFT-HR-R) 


0.137E+01 
0 . 137E+01 
0. 137E+01 
0.136E+01 
0.136E+01 
0. 136E + 01 
0.135E+01 
0. 134E+01 
0.134E+01 
0 . 134E+01 
0 .134E+01 
0 . 134E+01 
0 .134E+01 
0 . 1 34E+01 
0 . 134E+01 
0 . 134E+01 
0.134E+01 
0 . 1396 + 01 
0.134E+01 
0 . 134E+D1 
0.134E+01 
0 . 134E+01 
0 . 134E + 01 
0. 134E + 01 
0 . 134E+01 


0.Z54E+00 
0 . 255E + 00 
0 .255E+00 
0 .256E+00 
0 .257E+00 
0 .257E+00 
0 . 264E+00 
0. 268E + 00 
0. 269E+00 
0 . 269E+00 
0.269E+00 
0 .269E+00 
0 .270E+00 
0 .271E+00 
0 .272E + 00 
0 . 27 3E + 00 
0.273E+00 
0. 273E + 00 
0 . 269E + 00 
0 . Z70E+00 
0 .271E+00 
0 .Z71E + 00 
0 .272E+00 
0.271E+00 
0.271E+00 


0 . 2 1 7 6” 04 
0 . 22 0E- 04 
0.220E-04 
0.2Z4E-04 
0 .227E-04 
0.227E-04 
0 . 254E-04 
0 . 26 8 E- 04 
0 . 27 3E-04 
0 . 272E-04 
0 , 273E-04 
0 .273E-04 
0 .276E-04 
0 .279E-04 
0 . 282E-04 
0 .283E-04 
0.284E-Q4 
0 . 284 E- 04 
0 .272E-04 
0 .274E-Q4 
0 .276E-04 
0 .279E-04 
0 . 280E-04 
0 . 27 9 E - 04 
0.279E-04 


0 . 285E-01 
0 .289E-01 
0 . 291E-01 
0 .296E-01 
0 . 301E-01 
0 . 30LE-01 
0 . 348E-01 
0.373E-01 
0 . 382E-01 
0.381E-01 
0.382E-01 
0 . 382E-01 
0 . 387E-01 
0 . 394E-01 
0. 399E-01 
0 , 401E-01 
0 . 403E-01 
0 . 403E-01 
0 . 381E-01 
0 . 389E-G1 
0 . 388E-01 
0 . 393E-0L 
0. 396E-01 
0 . 393E-01 
0 .394E-01 


I RE 


PR ST 


MU 


F 


1 

0. 1640E+05 

0 .6959E+0D 

2 

0. 1118E+05 

0 .6960E+00 

3 

0 . 1259E+05 

0 .6961E + 00 

4 

0. 1465E+05 

0.695ZE+00 

5 

0 . 1465E+05 

0 . 6952E + 0G 

6 

0 . 1803E+05 

0 . 6926 E + 0 0 

7 

0 .2243E+05 

0.6932E+00 

8 

0 .2202E+05 

0 .6931E+00 

9 

0 .2937E+05 

0 .6931E+00 

10 

0 . 226ZE+05 

0 .6931E+00 

11 

0 . 2007E+05 

0 . 6931 E+ 0 0 

12 

0 . 1985E+05 

0 .6931E+00 

13 

0 . 1829E+05 

0 .6933E+00 

14 

0. 1258E+05 

0 . 69 36 E+ 0 0 

15 

0. 1013E+05 

0.6938E+OO 

16 

D . 8478E+04 

Q.6939E+Q0 

17 

0.8476E+04 

0 .6939E+00 

18 

0 . 1 029E+05 

0 .6931E+00 

19 

0 . 1023E+05 

0.6931 E+0Q 

20 

0 . 94456+ 04 

D . 6 932E + 00 

21 

0 -9163E+04 

0 .6933E+00 

22 

0 .8279E+04 

0 .6934E+00 

23 

0 . 1664E+05 

0 .6933E+00 

24 

0 . 1411E+05 

0 .6933E+00 


0 .4178E-02 
0 .4366E-D2 
0 4473E-02 
0 .4347E-02 
O.OOOOE+OD 
0 . 1607E-01 
0.1507E-01 
0. 1486E-01 
0 . 3853E-02 
0. 2898E-02 
0 . 3998E-02 
0 . 1532E-01 
0 . 1560E-01 
0 . 1706E-D1 
0 . 1832E-01 
0 .1934E-01 
0 . OOOOE+OO 
0 . 1168E-0I 
0 . 1170E-01 
0 . 1174E-01 
D.1177E-01 
0.1179E-01 
0 . OOOOE+OO 
0 . 4553E-02 


0 . 4978E+02 
0 4171E+0Z 
0 .3784E*02 
0 .4243E + 02 
0 . OOOOE + OO 
0 1876E+03 
0 .2202E+03 
0 2279E+03 
0 .6864E+02 
0.5224E+02 
0. 5918E + 02 
0.2114E+03 
0.2019E+03 
0. I619E+03 
0 . 1 356 E + 0 3 
0 . 1188E + 03 
D 00D0E+00 
0 .9119E + C2 
0 .9081E + 02 
0 . 9021 E+ 02 
0 . 896ZE + 02 
0 .89336+02 
0 . 0000E+00 
G . 4656E + 02 


0 .7204E-02 
0 .7565E-02 
0.7742E-02 
0 . 7532E-02 
0 . OOOOE + OO 
0 . 3043E+00 
0 . 3043E+0 0 
0 . 3043E + 00 
0 15I7E+OI 
0 . 1104E+00 
0 . 1467E+01 
0 . 3043E+ 00 
0 . 3043E+00 
0 . 3043E+00 
0 . 3043E+00 
0 .3043E+00 
0. OOOOE+OO 
0.9000E-01 
0.9000E-01 
0. 9000E-01 
0.9000E-01 
0 .9000E-01 
0 . OOOOE+OO 
0 .6991E-02 


FCOR 

HC 

HCCDR 

( BTU/FT*#2-HR-F) 


0. 1000E+01 

0. 9022E+02 

1 . 0000 

0 . 1000E+01 

0 .8976E+02 

1 .0000 

0 . 1OO0E+OI 

0 . 1020E+03 

1.0000 

0 . 1 000E+01 

0.1174E+03 

1 .0000 

0 . 1 000E+01 

0 . OOOOE+OO 

1 . 0000 

0 . 1 000E+01 

0 .9641E+03 

1 . 0000 

0 . 1000E + 01 

0 . 1300E+04 

1 . 0000 

0 . 10 0 0 E+Q 1 

0.13S6E+04 

1 .0000 

0. 1 000E + 01 

0 .4410E+03 

1 . 0000 

0 . 1 000E+01 

0 . 3547 E+ 03 

1 . 0000 

0 . 1000E + 01 

0 . 3932E+03 

1 . 0000 

0. 1000E + 01 

0 . 1383E+04 

1 . 0000 

0 . 1000E + D1 

0 . 1 305E+04 

1 . 0000 

0 . 1000E+01 

0. 9703E+03 

1 . 0000 

0 . 1000E+01 

0.7729E+03 

1 . 0000 

0 . 1000E+01 

0.6588E+03 

1.0000 

0 . 1000E+01 

0 .OOOOE+OO 

1 .0000 

0 . 1000E+01 

0. 3209E+03 

1.0000 

0 . 1000E+01 

0.316 1 E+03 

1.0000 

0. 1 0O0E+01 

0 .31596+03 

1 .0000 

0. 1000E+01 

0 .3157E+03 

1.0000 

0. 1000E+01 

0 . 3162E+03 

1.0000 

0. 1000E+01 

0 . OOOOE + OO 

1 . 0000 

0. 1000E+01 

0.257 5E+ 03 

1 . 0000 


HCCORC 

HTAC0R 

TAW 


(F) 

1.0000 

1.0000 

0 .76466+03 

1 .0000 

1 . 0000 

0.7734E+03 

1 .0000 

1 . 0000 

O.7989E+03 

1 .0000 

1 . 0000 

0 . 8245E+03 

1 . 0000 

1 , 0000 

0 . 0000E+00 

1.0000 

1 . 0000 

0 . 10Z5E+04 

1.0000 

1 . 0000 

0 . 1 186E+Q4 

1.0000 

1.0000 

0 . 1255E+04 

1 . 0000 

1.0000 

0 .1267E+04 

1 . 0000 

1 .0000 

0 . 1268E+04 

1.0000 

1.0000 

0 .1269E+04 

1.0000 

1.0000 

0 . 1294E+D4 

1.0000 

1.00 0 0 

0 . 133GE+04 

1 . 0000 

1.0000 

0 . 1360E+04 

1 .0000 

1.0000 

0.1374E+04 

1.0000 

1.0000 

0 .1385E+04 

1 .0000 

1 . 0000 

0 . 0000E+00 

1.0000 

1.0000 

0 . 1324E+04 

1.0000 

1 . 0000 

0 . 1275E+04 

1 . 0000 

1 . 0000 

0 . 1302E+04 

1 .0000 

1 . 0000 

0.1S31E+04 

1.0000 

1 . 0000 

0 . 1357E+04 

1.0000 

1 . 0000 

0 . OOOOE+OO 

1.0000 

1 . 0000 

0 . 1366E+04 



BYPASS FLOW BETWEEN INTERVALS 

OR SLICES 

2 AND 18 







(CONVERGED IN 

3 ITERATIONS) 

LBM/S 








FLOW RATE 3 

0 . 16056E-02 








DX = 

0 .65000E+00 

IN 








DHYD = 

0 .75000E-01 

IN 








FR. FCT = 

0 . 10000E+00 









IHLT. PR. REC . = 

0 . 50000E+00 









A 

( IN**2) 

R 

( IN) 

TT 

(F) 

T 

(F) 

PT 

(PSIA) 

P 

(PSIA) 

M 

V 

(FT/S) 

ALPHA 

(DEG) 

IN 

OUT 

0 . 1000E-02 
0 . 1000E-02 

0 . 3800E+01 0 

0.3750E+01 0 

.7733E+03 
. 8 099E+ 03 

0 . 7413E+05 
0 . 7 157 E+ 0 3 

0 . 1775E+03 
0. 1223E+03 

0 . 1609E+03 
0 .91756+02 

0 . 3809E+00 
0 .6610E+00 

0.6395E+03 

0.1098E+04 

0 .9000E+D2 
0 .3600E+02 


Figure 9.— Concluded. 


Note that a value of Z, R, Y, A, DH, and KT is shown before 
the numbered node or interval values. These are the inlet 
values of position, area, hydraulic diameter, and KTIN, 
respectively. 

A complete description of the specified tip-cap geometry 


and the calculated tip-cap flow follows. The output is self- 
explanatory, but the user should be aware that the tip-cap flow 
calculations do not account for radius change (pumping) in 
the impingement or exit holes. This does not introduce any 
significant errors, since for any realistic design the 
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impingement and exit hole lengths will be kept as short as 
possible. The tip-cap How calculations also have no provisions 
for calculating reverse flow. If too high a backpressure is 
specified at the tip cap, the tip-cap flow will not converge. 
This will be noted by appropriate error messages up to and 
including the final overall flow iteration. The obtained flow 
solution, although meaningless for the tip-cap geometry, is 
the same as the solution that would have been obtained had 
no tip-cap flow been specified. 

The next tabulation presents calculated flow properties and 
physical properties. As before, the inlet values precede the 
numbered node values (all values in this table are node values). 
The tabulated parameters are 
I node number 

T static temperature 

TT total temperature 

P static pressure 

PT total pressure 

W coolant flow rate 

V velocity 

M Mach number 

GAMMA ratio of specific heats 

CP specific heat at constant pressure 

VISC. viscosity 

TH. CNDT. thermal conductivity 

Note that for rotating turbomachinery (RPM ?£ 0) the shown 
total temperature and total pressure, as well as velocity and 
Mach number, represent values in the rotating reference frame. 

The next tabulation is a mixture of flow and heat transfer 
parameters, along with specified correction factors. The 
tabulated parameters are 
I node or interval number 

RE node Reynolds number based on hydraulic 

diameter (DH) 

node Prandtl number 

interval Stanton number based on uncorrected HC 
interval Nusselt number based on uncorrected HC 


PR 

ST 

NU 

F 


specified or calculated interval Fanning friction 
factor 

FCOR specified interval friction factor correction 
HC specified or calculated interval heat transfer 

coefficient w ithout any correction factors 
HCCOR specified interval heat transfer coefficient cor- 
rection factor 

HCCORC calculated interval heat transfer correction factor 
due to rotation (listed as 1 .0 if IOPTCR(/ 2 > = 0) 
HTACOR specified interval heat transfer area correction 
factor 

TAW calculated interval adiabatic wall temperature 
(listed as 0.0 for fake intervals) 


Further explanation is in order for the following calculated 
interval output parameters: ST, NU, F, HC, HCCORC, and 
TAW (excluding F and HC if they are user specified). If the 
slice option is used (NSLI(rt) > 1) together with the option 
for printing output at all the slices (ISLOUT = 1), it is obvious 
that the shown interval values correspond exactly with the slice 
interval numbers. However, if the slice option is used together 
with the option for printing output in the original frame of 
reference (ISLOUT = 0), an arbitrary decision must be made 
as to which of the calculated slice interval values (the slice 
intervals that were obtained from the original input interval) 
will be printed. For the described combination of options 
(NSLI(w) > 1 and ISLOUT = 0), the shown calculated 
interval values (ST, NU, F, HC, HCCORC, and TAW) are 
the values for the last (highest numbered) slice interval in the 
original input interval. 

The final output is a tabulation of bypass flow parameters 
including a listing of values at the bypass inlet and exit. This 
output is self-explanatory, except for the following item: CPF 
has no provision for calculating reverse flow in the bypass 
channel. If the main coolant flow cannot be balanced without 
reverse flow in the bypass, the bypass flow is set to zero and 
the main coolant flow calculations are carried through to 
convergence. The bypass flow output will show zero flow, 
as well as zero velocity and Mach number. The obtained 
solution, although meaningless for the specified bypass 
geometry, is the same as the solution that would have been 
obtained had no bypass been specified. The user is therefore 
cautioned to always check that the calculated bypass flow rate 
is greater than zero. 


Error Messages 

All DO loops in CPF have appropriate error messages that 
are printed if convergence is not achieved. A detailed listing 
of the error messages will not be provided here, since most 
messages are self-explanatory and contain a suggested course 
of action to alleviate the problem. The follow ing comments, 
however, should be helpful in understanding the error messages: 

(1) All error messages stating node or interval numbers 
refer to the calculated slice frame of reference. 

(2) Error messages pertaining to bypass or tip-cap flow are 
often caused by specifying too large a flow area. 

(3) Error messages fall into two general groups: warning 
messages and abort messages. Program execution continues 
after a warning message is issued and terminates after an abort 
message. 

(4) Although not all warning messages use the word 
“WARNING,” all messages that result in program termi- 
nation use the word “ABORT.” 

(5) Most warning messages have the following generalized 
format: 

WARNING— XXXXX HAS NOT CONVERGED IN 
OVERALL FLOW ITERATION XX 
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These messages can be triggered by extreme initial conditions 
during the iteration process. The solution is valid if the warning 
messages are for an overall flow iteration number XX that 
is less than the required number of total overall iterations 
NTOT. This number (NTOT) is shown following the listing 
of any warning messages in the following statement: 

CONVERGENCE ACHIEVED IN NTOT 
OVERALL ITERATIONS 

(6) The user is cautioned that it is possible to specify 
configurations or situations that have no physical solution. 
These are usually caused by specifying the flow option 
IOPTW = 1 while also specifying incompatible inlet and exit 
flow areas and pressures, as well as possibly unrealistic friction 
factors or heat transfer coefficients. For such cases the 
following error message is printed out: 

ABORT— COOLANT FLOW HAS NOT CONVERGED 
IN 100 ITERATIONS 

For this error message a suggested course of action is described 
under Miscellaneous Comments in the section Program Input. 


Example Problems 

The example problems shown in this section were run on 
an IBM-compatible Amdahl-5869 computer and with an initial 
version of CPF which had no provision for direct heat addition 
(see example 3). 

Example 1 

Example 1 illustrates most of the features of CPF by 
analyzing the upper flow path of the conceptual cooled radial 
turbine of Figure 1. This figure does not represent a specific 
cooled radial turbine design but illustrates the various features 
that CPF can accommodate. 

The nodal breakdown of the upper flow path (consisting of 
24 nodes) is shown in figure 2, and the detailed data input is 
shown in figure 8. The plain passage (from the inlet to the start 
of the turbulators) is modeled by intervals 1 to 4. Specified 
options for these intervals are IOPTHT = 4 (internal plain- 
passage correlation) and IOPTFF = 2 (user-input friction factor 
curve using NPFF4, etc.). Note that interval 2 contains the inlet 
for the bypass flow (NBPI = 2). Interval 5 is a fake interval 
that accounts for the sudden area change between the plain 
passage and the trip-strip passage. Dummy values (equal to 7) 
are used for IOPTHT and IOPTFF. Intervals 6 to 8 model the 
trip-strip passage up to the first turn at the tip cap (IOPTHT = 1 
and IOPTFF =1). Interval 9 models the 90° turn via a plain- 
passage heat transfer correlation (IOPTHT = 4) and a specified 
AT factor (KTSGMT) for friction (IOPTFF = 4). In interval 10 
the heat transfer is specified from an input Colburn /-factor 
curve (IOPTHT = 5 and using NPCJF, etc.), along with a 


specified K factor for friction (IOPTFF — 4). Interval 10 also 
feeds the tip-cap How (NTCF = 10). Interval 11 models the 
second 90° turn at the tip in the same manner as interval 9 
(IOPTHT = 4 and IOPTFF = 4). Intervals 12 to 16 model the 
trip-strip intervals leading up to the start of the pin tins 
(IOPTHT = 1 and IOPTFF =1). Interval 17 is a fake interval 
that models the sudden area change between the turbulators and 
the pin fins. As before, a dummy value of 7 is specified for 
IOPTHT and IOPTFF. Intervals 18 to 22 model the pin-fin 
passage up to the finned exit (IOPTHT = 2 for heat transfer, 
along with a user-specified friction factor, IOPTFF = 0). Note 
that the bypass flow reenters at interval 18 (NBPO = 18). 
Interval 23 is a fake interval to account for the sudden area 
change at the finned passage (a dummy value of 7 for IOPTHT 
and IOPTFF). Interval 24 models the finned exit passage (using 
IOPTHT = 3 and IOPTFF = 3). 

The following discussion, as well as the notes on the left 
margin of figure 8 further describe and clarify the input. The 
first line is the title and the second line starts the NAMELIST 
data set called DATT. The first three lines set overall param- 
eters as well as inlet and exit conditions. Note that English 
units are specified and that two slices are specified for each 
of the 24 intervals (NSLI = 24*2). The next four lines specify 
parameters for the input Fanning friction factor curve, and 
the following three lines specify the Colburn /-factor input 
curve. The geometric flow path (at the passage centroid) is 
described via the parameters Z, R, and Y. There are 24 
specified values for each, followed by 126 zeros. These trailing 
zeros are arbitrary values that further fill the arrays and have 
no bearing on the calculations. Such till-in values are used 
throughout the data input set. The next four lines are the 
specifications of the options for fake intervals, heat transfer, 
and friction factor. Note that, where dummy values are needed, 
the value of 7 is used. The next two lines contain geometric 
parameters (EDD and PDE) needed for the specified trip-strip 
heat transfer correlation (IOPTHT = 1), the roughness 
parameter (DOE) tor the specified friction factor correlation 
(IOPTFF = 3), and specified total pressure loss coefficients 
(KTSGMT). Note that although these values are needed only 
for specific intervals, values are assigned for the total array. 
The next four lines contain parameters needed for the specified 
finned passage and are followed by eight lines with required 
area and hydraulic diameter values. Note that the parameters 
AV and DHV are needed only for the pin-fin passages. The 
next line contains values of user-specitied friction factors (F) 
and wall temperatures (TW), and the last six lines specify the 
needed parameters for tip-cap and bypass flow. 

The CPF output corresponding to the figure 8 data input 
is shown in figure 9. The output is in the original node/interval 
input frame of reference. Note that warning messages appear 
for overall flow iterations 8 and 9. The warnings can be 
ignored, since convergence was achieved in 27 overall 
iterations. Also note that for fake intervals the output values 
of Stanton number (ST), Nusselt number (NU), friction factor 
(F), heat transfer coefficient (HC), and adiabatic wall 


ll 


temperature (TAW) are all shown as zero. Finally, note that, 
although each interval was specified to be divided into two 
slices (NSLI = 24*2 in fig. 8), figure 9 shows that for the 
appropriate intervals (2. 5, 9. 10, 11, 17, 18, and 23) no 
interval division was allowed (NSLI forced to be 1 for fake 
intervals, bypass intervals, tip-cap intervals, and intervals with 
specified KTSGMT). 

Example 2 

As stated in the section Flow Path Modeling, the free flow 
area change between nodes or slices defining a “regular” 
interval determines the accuracy of the flow calculations. 
Example 2 illustrates the amount of computational error to be 
expected from various area variations. Two flow paths are 
analyzed: the first decreasing in area from 10 cm 2 to 5 cm 2 , 
the second increasing in area from 5 cm 2 to 10 cm 2 . Both 
flow paths are modeled with just one interval, which is 


successively divided into various slice numbers (from 2 to 
150). The analysis assumes no friction and no heat transfer 
(isentropic flow) in order to isolate the effect of area change. 
Physical properties are assumed to be constant. The input data 
are shown in figure 10. Note that hydraulic diameter (DH) 
is set to a dummy value (1.0), since no friction or heat transfer 
is specified. Inlet total pressure and total temperature are 
identical for both cases (100 N/cm“ and 500 K, respectively); 
exit pressures are 65 and 94 N/cni“ for the decreasing and 
increasing area cases, respectively. The fixed weight flows 
correspond to the specified inlet total and exit static conditions. 
The calculated results are summarized in table II. 

Since no friction or heat transfer is specified, the exit total 
conditions ideally should be equal to the inlet total conditions. 
As can be seen, this is not the case, and the changes in total 
pressure and total temperature from inlet to exit are thus a 
measure of the flow calculation accuracy. As expected, 
accuracy increases with an increase in the number of specified 


TABLE It. -SUMMARY OF EXAMPLE 2 RESULTS (OUTLET CONDITIONS) 
(a) Decreasing flow area 


A in = 10 cm 2 
Pm - 100 N/cm 2 
7in - 500 K 





CPF predictions for number of slices 

Insentropic 

solution 

2 

5 

10 

20 

40 

60 

80 

100 

125 

150 

Temperature, T, K 
Total temperature, T\ K 
Pressure, p y N/cm 2 
Total pressure, /?', N/cm 2 
Mach number, M 

— 

439.1 

497.3 

63.69 

99.27 

0.8247 

441.7 

498.6 

64.81 

99.82 
0.8128 

442.5 

499.0 

65.13 

99.98 

0.8094 

442.7 

499.1 

65.23 

100.0 

0.8085 

442.7 

499.1 

65.24 

100.0 

0.8083 

442.7 

499.1 

65.25 

100.0 

0.8083 

442.7 

499.1 

65.25 

ioo.o 

0.8082 

442.7 

499.1 

65.25 

100.0 

0.8083 

442.7 

499.1 

65.25 

100.0 

0.8083 

443.1 

500.0 
65.0 

100.0 
0.8117 


(b) Increasing flow area 


A m = 5 cm 2 
p' n = 100 N/cnr 

r; n = 5oo k 



10 cm 2 


= 94 N/cnr 




CPF predictions for number of slices 

Insentropic 

solution 

2 

5 

10 

20 

40 

60 

80 

100 

125 

150 

Temperature, T, K 

500.5 

495.1 

493.3 

492.6 

492.4 

492.3 

492.3 

492.3 

492.3 

492.2 

491.4 

Total temperature, T , K 

508.9 

503.6 

501.9 

501.2 

501.0 

501.0 

501.0 

500.9 

500.9 

500.9 

500.0 

Pressure, p y N/cnr 

97.39 

95.12 

94.36 

94.06 

93.96 

93.94 

93.94 

93.93 

93.93 

93.92 

94.0 

Total pressure, />', N/cnr 

103.3 

101. 1 

100.4 

100.1 

99.98 

99.96 

99.95 

99.95 

99.94 

99.94 

100.0 

Mach number, M 

0.2919 

0.2973 

0.2991 

0.2999 

0.3001 

0.3001 

0.3002 

0.3002 

0.3002 

0.3002 

0.2997 
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EXAMPLE PROBLEM 2 (CALCULATION ACCURACY) 
g DATT IECHO=0, IUNITS=1, IOPPRP=0, RG=287.053, RPM=0.0, 

G= 1 . 39 , CP=0.24, XMU=2.65E-4, XK=0.95E-4, 

NI=1,ZIN=0 . 0, RIN = 0.0» Y I N = 0 . 0 , AIN=10.0, DHIN-l.O, 

PTSP = 1 00 . , TTSP = 500 . , KTIN=0.0, POUT=65.0, 

Z-10.0, R=0.D, Y= 0 . 0 , A = 5 . 0 , IOPTW=0, NIN=0.8721, 

I0PTFI = 0, I OPTHT = 0 , IOPTFF=0, N$LI=2, 

HC = 0 . 0 > DH=1.0, F=0.0, TN=1000. SEND 

&DATT NSL 1=5 SEND 

SDATT N S L I = 1 0 SEND 

SDATT N S L I = 2 0 SEND 

SDATT NSL I =40 SEND 

SDATT NSL I =6 0 SEND 

SDATT N S L I = 8 0 SEND 

SDATT NSL I = 1 0 0 SEND 

SDATT NSL I = 125 SEND 

«DATT AIN^slo? POUT=9A.O, A=10.0, WIN=0.88A3, N$LI=2 SEND 

SDATT NSL I =5 SEND 

SDATT NSL I =1 0 SEND 

SDATT NSL I =20 SEND 

SDATT NSL I =40 SEND 

SDATT NSLI=60 SEND 

SDATT NSL I =80 SEND 

SDATT NSL I = 1 0 0 SEND 

SDATT NSLI =1 25 SEND 

SDATT NSLI=150 SEND 

Figure 10.— Example 2 input data. 


slices. A two-to-one area change from inlet to exit is a severe 
gradient, and a small number of slices cannot yield an accurate 
solution. Note that for the decreasing-area case, no solution 
was obtained for two slices (the calculated exit Mach number 
exceeded 1.0). Table II shows that, for both cases, solution 
accuracy increased up to approximately 40 slices but remained 
essentially unchanged for slice numbers greater than 40. 

The fact that the CPF predictions do not match the exact 
solution is not surprising, since CPF is by nature ill suited 
for calculating isentropic flow paths (owing to the extensive 
nested iteration loops in the calculation procedure). The results 
suggest that sufficient engineering accuracy is obtained if the 
free flow area change between adjacent nodes or slices is less 
than 10 percent. 

Example 3 

Example 3 compares the CPF predictions for both Fanno 
line and Rayleigh line calculations for Mach numbers between 
0.20 and 0.75. Both cases treat a constant-area duct with an 
inlet area of 10 cm 2 (diameter, 3.568 cm) and an inlet total 
temperature and pressure of 300 K and 100 N/crrr, 
respectively. Weight flow is fixed at 0.7874 kg/s, which results 
in an inlet Mach number of 0.20. The flow paths are modeled 
with one interval and various numbers of slices (50, 100, and 
150). The input data are shown in figure 11. 

A friction factor of 0. 10 along with zero heat transfer was 
chosen to check the Fanno line prediction. From standard 
Fanno line tables the needed path length to go from Mach 0.20 
to 0.75 (with a friction factor of 0.10) is 36.0143 diameters 
(Z = 36.0143*3.568 = 128.50803 cm in fig. 1 1). Calculated 
results are shown in table III for specified slice numbers of 
50, 100, and 150. Calculation accuracy increased with slice 
number, and for 150 slices the predicted results essentially 
coincided with the Fanno line predictions. 


EXAMPLE PROBLEM 3 (FANNO AND RAYLEIGH LINE) 

SDATT IECH0=0, IUNITS=1, I0PPRP=0, RG=287.053, RPM=0.0, 
G= 1 . AO , CP=0. 23978, XMU=1.8E-4, XK=0.6E-4, 

NI = I, Z I N - 0 . 0 » R I N = 0 . 0 , Y I N = 0 . 0 > AIN = 10.0, DHIN = 3 . 56825, 
I0PTW=0, HIN=. 787355, 

PTSP= 100.0, TTSP=300 . 0, KTIN=0.0, 

Z=128. 50803, R-0 . 0 , Y=0.0, NSLI=50, A=10. 

F = 0 . 10 , TW=500 . SEND 
SDATT NSL I =1 00 SEND 

SDATT NSLI = 150 SEND oeiin 

SDATT Z = 12 . 576 , F=0.0, HC=2.0, TW=2000., NSLI=50 SEND 
SDATT Z = 1 2 . 57 A , NSLI = 100 SEND ' 

SDATT Z = 12 . 57 2 , NSLI = 15G SEND 


HC = 0 , Q , DH = 3. 56825, 


Figure 11. — Example 3 inpul data. 


TABLH III. -SUMMARY OF EXAMPLE 3 RESULTS 


(a) CPF Fanno line predictions. Friction factor. 0.10; path length. 
36.0143 diameters. 



Inlet 

conditions 


Outlet 

conditions 

CPF prod id ions for 
number of slices 

Fanno 

solution 

50 

100 

150 

T emperature, T, K 

297.6 ! 

268.9 

269.5 

269.6 

269.6 

Total temperature, 7’\ K 

300.0 

299.7 

299.9 

299.9 

300.0 

Pressure, />, N/cnr 

97.25 

24,46 

24.66 

24. 68 

24.69 | 

Total pressure, p\ N/cur 

100.0 

35.72 

35.83 

35.84 

35 . 85 j 

Mach number, M 

0.20 

0.756 

0.751 

0.750 

0.750 | 


(h) CPF Rayleigh line predictions. Arbitrary heal transfer coefficient and uall 
temperature; path length varied to produce desired exit Mach number. 



Inlet 

conditions 

Outlet conditions j 

CPF predictions for 
number of slices 

Rayleigh 

solution 

50 

100 

150 

Temperature, 7*. K 

297.6 

1461 

1461 

1460 

1461 

Total temperature, T ' , K 

300.0 

1625 

1625 

1625 

1625 

Pressure, p, N/cirr 

97.25 

57.45 

57.45 

57.45 

57.45 

Total pressure, p' , N/eirr 

100.0 

83.43 

83.43 

83.43 

83.44 

Mach number, M 

0.20 

0.750 

0.750 

0.750 

0.750 

Path length, cm 


12.576 

12.574 

12.572 

— 


Closed-form solution Rayleigh line results arc based on fixed 
levels of heat flux into the fluid. Since the initial version of 
CPF (by which this example was calculated) did not accommo- 
date specified levels of heat flux directly, comparison of results 
had to be based on factors other than total heat flux. For the 
CPF calculations an arbitrary heat transfer coefficient and wall 
temperature were specified (HC = 2.0 cal/cirr s K and 
TW = 2000.0 K), along with zero friction factor. Path length 
was varied until the desired exit Mach number (0.75) was 
reached. The input data are shown in ligure 1 1 (last three 
lines). A comparison of the calculated exit temperatures and 
pressures (for the various slice numbers) is shown in table III. 
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All three cases match the Rayleigh solution temperatures and 
pressures. Note, however, that slightly different path lengths 
had to be specified for the various slice numbers. Appendix E 
describes a feature added to CPF by which the Rayleigh line 
problem can be solved directly. A heat generation term has 
been added to the energy equation; an example problem is 
illustrated for both positive and negative heat generation, and 
agreement is shown with the closed-form solution. 

Concluding Remarks 

A computer program called Coolant Passage Flow (CPF) 
has been developed that calculates the coolant’s flow rate, 
temperature, pressure, and velocity and the heat transfer 
coefficients inside turbomachinery cooling passages. Flow can 
be bled off for tip-cap impingement cooling, and a flow bypass 
can be specified in which coolant flow is taken off at one point 
in the flow channel and reintroduced at a point farther 
downstream in the same channel. Although specifically 
developed for radial turbomachinery, CPF can be used to 


analyze any turbomachinery coolant geometry that consists of 
a single flow passage with a single inlet and exit. CPF 
integrates the one-dimensional momentum and energy 
equations along a defined passage, taking into account area 
change, mass addition and subtraction, pumping, friction, and 
heat transfer. CPF predictions closely agree with closed-form 
solutions for isentropic flow, Fanno line flow, and Rayleigh 
line flow. Although the present version of CPF does not allow 
generalized flow branching, this ability is being developed. 
Reference 5 is a preliminary report describing an extension 
of CPF that includes a larger number of available heat transfer 
and friction factor specification options. 

For the flow passage geometries that it is designed to 
accommodate, CPF will give quick and accurate solutions. 
It is hoped that turbomachinery designers will find CPF to be 
a useful tool. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 26, 1990 
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Appendix A 



Symbols 




V 

velocity, m/s; ft/s (or open volume in 

A 

area, m 2 ; ft 2 


pin-fin interval, m 3 ; ft 3 ) 

Q> 

discharge coefficient; defined as ratio of actual 

w 

coolant flow rate, kg/s; Ibm/s 


to ideal flow 

X 

distance along flow path, m; ft 

Cl, C2, CT 

constants 

y 

distance in plane perpendicular to both r 

C p 

specific heat at constant pressure. 


and z (fig. 1). m; ft 


cal/(g K); Btu/ (Ibm °R) 

z 

distance along axis of rotation (fig. 1), m; ft 

D 

diameter, m; ft 

a 

angle between flow path and a perpendicular 

D h 

hydraulic diameter, 4/1/P, m; ft 


line to axis of rotation, deg 

e 

surface roughness, m; ft 

7 

ratio of specific heats 

e + 

parameter defined by e + = Re(W D) \ff/ 2 

A 

recovery factor; equal to V Pr or (Pr)° 333 

f 

8c 

Fanning friction factor, ApD h g c /2(V~p 
force-mass conversion constant, 1 ; 

M 

(eqs. (B29) and (B30)) 
viscosity, kg/(m s); lbm/(ft s) 
density, kg/m 3 ; lbm/ft 3 

32. 174(lbm ft)/(lbf s 2 ) 

P 

h 

enthalpy, J/kg; Btu/lbm 

T f 

fin thickness (fig. 6), m; ft 

/ - — \ o .5() 

h. 

coolant heat transfer coefficient, 

<Pf 

parameter defined by = \ 2h c ?? / kr/j 

J 

cal/(m 2 hr K); Btu/(ft 2 hr °R) 
Colburn 7-factor, h c (Pr) m /pVC p 


(eq. (B45)) 

K T 

total pressure loss coefficient ( K factor). 

Subscripts: 



Ap' f(U2pV 2 ) 

aw 

adiabatic wall 

k 

thermal conductivity, cal/(m hr K); 

b P 

bypass 


Btu/(ft 2 hr °R) 

ch 

choke 

L 

total length of flow path to a particular 

F 

film 


node, m; ft 

imp 

impingement 

h 

half-height of a fin (fig. 6), m; ft 

in 

inlet 

t 

interval length, m; ft 

tnj 

injection 

M 

Mach number 

n 

node or interval number 

m 

fin spacing (fig. 6), m; tt 

PP 

plain passage 

N 

rotational speed. 1/s 

ref 

reference 

Nu 

Nusselt number, h c /D h k = (St)(Re)(Pr) 

S 

supply 

P 

perimeter, m (ft) 

tc 

tip cap 

Pr 

Prandtl number; C p p/k 

V 

based on volume 

P 

pressure, N/m 2 ; psia (or trip-strip spacing 

w 

wall 


(fig. 5), m; ft) 

X 

component in direction of main-stream flow 

R 

gas constant, J/(kg K); (ft lbf)/(lbm °R) 

2 

station at tip-cap inlet 

Ra 

Rayleigh number, rN 2 (1 IT)(T W — T. dW )Dl 
X Pr /(p/p) 2 

3 

station in tip-cap chamber 

Re 

Rs 

Reynolds number, pVD b /p 
Rossby number, ND h /V 

4 

station at tip-cap exit 

r 

radius, m; ft 

Superscripts: 


S 

total heat transfer area of pin-fin interval. 

O' 

total conditions 


m ; ft 

( ) 

arithmetic average 

St 

Stanton number, Nu/(Re)(Pr) — h c /pVC p 

/ 

interval 

T 

temperature, K; °R 

nd 

nodal 
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Appendix B 

Analysis Equations and Flow Solution Procedure 


All flow calculations in CPF are based on one-dimensional, 
compressible fluid flow for a gas obeying the perfect-gas law. 

The following notation will be used throughout appendix B 
to differentiate between values at node points and interval 
values: symbols representing nodal values will be shown as 
plain symbols without modifiers: if the symbol refers to a 
specific node number, the symbol will be subscripted by the 
node number. For example. A, T, and C p represent general 
nodal values of area, temperature, and specific heat ratio, 
respectively; A n , /,,, and C p „ designate those parameters for 
node n. All symbols representing interval values, including 
those designating the arithmetic average of bounding nodal 
values, will be shown with a bar (~) over the symbol (denoting 
average value over the interval). Values for specific intervals 
will again be subscripted by the interval number. 

Main Program 

The main program starts by setting default values and by 
checking for incompatibilities between the selected friction and 
heat transfer options. The specifications for tip-cap and bypass 
flow are also checked. If any irregularities are found, the 
program terminates with an appropriate error message. 

The chosen input units are next converted to working English 
units, and the maximum allowable (choked) weight flow at 
the inlet is calculated (to establish an upper bound on allowable 
weight flow) as follows: 


Th = ■ 


1 + 


T[ 

7 - 


Pch — 


Kk = (yftg,T th )°- 

Ps 


50 


1 + 


7 - 1 


yf(y- M 


Pc h 


Ps - ^’r.infl.-h^ch 

2RgJ h 

Pch K-'h^in 
R gJ , ’ ch 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 


flow path from node to node. If the user chooses to subdivide 
the input flow intervals into slices (NSLI(n) > 1 ), the program 
switches from the input nodal frame of reference to a new, 
calculated, slice frame of reference, in which new node and 
interval numbers are created in accordance with the total 
number of slices. All calculations are performed in the new 
frame of reference. The required node parameters needed at 
the new slice node locations are obtained by linear interpolation 
of the appropriate original input nodal values. Input variables 
that were specified for the original input intervals are retained 
for the new intervals (which were created by slicing up the 
original intervals). For all equations presented in this appendix, 
references to node or interval numbers can refer to either the 
original input frame of reference or the new, calculated, slice 
frame of reference, as appropriate. 

The pressure and temperature at each node are calculated 
from the previous node by 


Pn 


Pn — I + 



(B 6 ) 


T n = Tn i + (~^j dx n (B7) 

The derivatives (■ dp!dx) n and (dT!dx) n are derived in appendix D 
and are as follows: 



+ 




V wCp 2wC p ) 


+ Vinj* _ 2(w/A)Rg c T ") /dw\ 

A p a j \*/ B 


((w/A) 2 Rg c T 


pA 


Cl / wRg c T \ 2 '} / dA \ 
C p A \ pA ) 1 \dx) 

i N / J \ / n 


If the physical properties are calculated from subroutine 
AIRPRP , these calculations are iterative, since 7 is then a 
function of temperature. 

Overall flow iteration procedure.— From the specified 
entrance conditions and an assumed or specified inlet flow rate, 
the flow solution is obtained by marching along the defined 


4fRg,-T(w/A) 2 N 2 rp / dr\ 
2pD h + RgJ \dx) 


4h, AC\(T w - 7~ aw ) 
wCpD h 


(B 8 ) 
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<rr\ _j_ 
dx) ~ ct 

/ n 


w (RgJX + 1 + i Rg ‘ T 


C n \ pA / vv 2 C p \ pA 


r T 

'-'/j.mj 1 mj 

C p w 


V\ 


,n J 


2 wC, 


dw\ 
dx ) 

/ n 


+ 


1 f wRgj X} ( dA\ 


C P A\ P A ) \dxj n 


+ 


1 f\vRg l fY\(dp 


dx 


1 pC p \ pA 


+ 


Ah,A( 7j, - 7; w ) 
wCpD h 


(B9) 


where K inj _ t is the component of velocity of the injected gas 
in the main-stream flow direction and 


dw\ _ w„ - w„-| 


dx ) n 


(BIO) 


dA \ A n A n ~ | 

dx) n („ 


(Bll) 


r , i ~ r n - 1 


( B 1 2 ) 


The constants Cl, C2, and CT are defined as follows: 


Cl = 


(w/AyRg c 


1 + 


(w/A) 2 R 2 g*T 
P 2 c P 


(B13) 


(w/A) 2 Rg c T Cl 
C2 = 1 — ' — 4- 


P c p 


(w/A)Rg t T 


(B14) 


CT= 1 + 


(w/A)*R z g;.T 

P 2 C P 


(B15) 


Equations (B6) to (B15) contain variables denoting average 
values over the interval from node n - 1 to node n. They 
are designated by a bar ( ) over the symbol. The average 
geometric parameters (r, A, and D/ z ) are known, but average 
flow and heat transfer parameters (vv, p , T, h, and T aw ) and 


average physical properties (C p ) must be determined in an 
iterative manner during the solution marching process. The 
formulation of equations (B6) to (B15) also places a premium 
on small step size for calculation accuracy. Although the slice 
option (NSL1 (aO > 1) can be used to generate any number of 
additional slices, each new node and interval generated by the 
slice option creates additional array requirements for numerous 
variables. Depending simply on sufficient slices to get the 
desired calculation accuracy requires very large array sizes 
for the program variables. 

In order to get additional accuracy without needing an 
exorbitant number of slices, the program further divides each 
input or slice interval into smaller sections called steps. New, 
average, geometric variables (position, area, and hydraulic 
diameter) are calculated by linear interpolation for each step. 
Pressure and temperature at the next node (node n) are 
determined from the previous node (node n - 1) values by 
the following iterative procedure, consisting of an inner and 
an outer iteration. 

For the first pass through the equations, all average values 
of pressure p and temperature T are assumed to be equal to 
the starting node (node n — 1) values. The interval is then 
divided into four steps and the appropriate average step 
geometric variables are calculated. The solution is marched 
to node w, where Mach number, pressure, and temperature 
are determined. Next, the number of steps is doubled, average 
step geometric variables are recalculated, and the solution is 
again marched from node n — 1 to node n. The calculated 
Mach number at node n is now compared w ith the previously 
calculated Mach number. The doubling ot step numbers is 
repeated until the relative difference between successively 
calculated Mach numbers at node n is less than 5.0 x 10 s 
(relative difference is defined as the absolute value of (New 
value— Old value)/01d value). This iteration procedure is 
referred to as the “inner iteration.’ ’ Note that neither physical 
properties, friction factors, nor heat transfer coefficients are 
recalculated at each step. They are based on average values 
in the interval. 

Once the Mach number has converged, the calculated 
pressure and temperature at node n are used to determine new 
average values of p and T between the nodes. The described 
inner iteration is repeated (for 4, 8, 16, . . ., etc., number 
of steps) with the new p and T values. This overall procedure 
is referred to as the “outer iteration. It is repeated until the 
relative differences between successively calculated values of 
p and T are less than 5.0X10' 4 . 

The described marching procedure is repeated for each 
interval from the inlet to the exit. If the flow rate was specified 
(IOPTW = 0) and all node Mach numbers (except at the exit) 
are less than 1 .0, the solution is obtained in one pass. If a Mach 
number > 1.0 is calculated, the program aborts with an 
appropriate error message. If the program is to determine the 
flow rate (IOPTW = 1) and a nodal Mach number > 1.0 
is calculated anywhere except at the exit, the marching 
procedure is restarted at the inlet with a reduced weight flow. 
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Once the exit is reached, the calculated pressure at the last 
node is compared with the specified exit pressure. A new guess 
of mlet flow rate is then made, and the marching procedure 
is repeated until the calculated and specified exit pressures 
converge. For choked flow at the exit the specified back- 
pressure is decoupled from the solution. 

When a converged solution is reached, the working English 
units are converted back to the specified input units and the 
output is printed. The input data set is then scanned for the 
next input case. 

Inlet conditions for bypass and tip-cap flows.— At the inlet 
to the bypass channel (in interval NBPI) the total pressure and 
temperature are a function of the assumed inlet pressure 
recovery factor bp (input variable PBPRI) and the average 
interval static and total pressures and temperatures as follows: 


/ 7 hp.m Pmm + (p NBPI / ; NBPl)^7.bp 


T — T 
1 bp. in 1 N 


+ ( 7 ' 


NBPI 


- 71 


NBPI 


)*tj 


bp 


where 


„ ' _ /' NBPI I - /'NBPI 

/'NBPI 1 


(B 16) 
(B 17) 


(B 18) 


„ _ / ; NBPI I + /'NBPI 

/'NBPI ~ 


(B 19) 


defined by 


An = 4/ 


/_r \ _^ 

\Dh) 2ft. 


(B24) 


Only five correlations for specifying friction or pressure drop 
in an interval are provided in CPF: (1) user-supplied friction 
factor , (2) trip-strip correlation, (3) friction factor input curves 
(tor pin fins, tinned passages, or plain passages), (4) roughened- 
p la in-passage correlation, and (5) /f-factor specification. No 
attempt is made to accommodate the myriad of available 
correlations, most of which have a narrow range of geometric 
applicability. The user can readily program specific correl- 
ations into CPF. Note that a broad range of correlations is 
used and discussed in reference 5. 

When user supplied, the friction factor is input via the input 
variable F(NI). When friction factor curves are input in tabular 
form as a function of Reynolds number for pin fins, finned 
passages, and plain passages, the corresponding input friction 
factor variables are FFF2(NPFF2), FFF3(NPFF3), and 
FFF4(NPFF4), respectively. For this option the required 
friction factor is determined through interpolation by sub- 
routine SPLFIT. which will be described separately. 

For trip strips (roughened tubes) the friction factor is 
calculated from a correlation of reference 6. 


7Vm>i — 


7 NBPI - I 4- 7 NHpj 


Tnbpi ~ 


7’nbpi i + 7^ B pj 


(B20) 

(B21) 


/ = 


2.0 


2.5 In 



- 3.75 + 0.95 


(B25) 


This equation is valid for 10 < p/e < 40 and e + > 35, where 


B Kj bp is zero, the bypass inlet total pressure and 
temperature are equal to the average static values in interval 
NBPI; it K } bp equals 1 .0, the bypass inlet total pressure and 
temperature are equal to the average total conditions in interval 
NBPI. The bypass inlet pressure recovery factor (default value, 
0.0) must assume a value between 0.0 and 1 .0. If a value less 
than 0.0 or greater than 1.0 is input, the variable is reset to 
0.0 and 1.0, respectively. 

The inlet total temperature and pressure for tip-cap flow in 
interval NTCF are assumed to be the arithmetic averages of 
the static conditions in the interval. That is. 


r 7’ _ Tntcf _ j 4- /n i c 

1 a . in “ * NTCF 

? 


_ /' NTCF -I + /'NTCF 
Ph in /'NTCF ~ “ 


(B22) 

(B23) 


Friction factor and pressure drop calculations .— The 
friction factor used in CPF is the Fanning friction factor 


e 


+ 




(B26) 


If the friction factor is to be calculated from the plain-passage 
correlations of reference 7, subroutine WWWFF is used. This 
subroutine will be described separately. 

If a total pressure loss coefficient K r is" specified for an 
interval (input variable KTSGMT), it is transformed to a 
Fanning friction factor as follows: 


f = 


KjD h 

4t 


(B27) 


Total pressure losses due to sudden expansions and con- 
tractions in fake intervals are calculated by subroutine 
TPLOSS, which will be described separately. 

Heat transfer calculations .— Heat addition to and subtrac- 
tion from the fluid in each interval are driven by the difference 
between the specified wall temperature and the calculated 
adiabatic wall temperature, which is defined by 


is 



r aw = t + A(r - d 


(B28) coefficient is then evaluated from 


The recovery factor A is a function of the average Reynolds 
number in the interval 

A = (Pr)° 50 for Re < 2300 (B29) 

A - (Pr) 0 - 333 for Re > 2300 (B30) 

Physical properties in the heat transfer calculations are 
evaluated at the reference temperature given by 


J, = (axiteXM ) ,836) 

D„ 

For pin fins (IOPTHT - 2) the correlation of Van Fossen 
(ref. 8) is used with hydraulic diameter and Reynolds number 
based on open volume. 

(Nu) v = 0. 153 (Re), 0 685 (B37) 


f Kf = 0.5(7’,,. — T) + T for Re < 2300 (B31) 

T nf = 0.5 7„. + 0.28 f + 0.22 T avi for Re > 2300 (B32) 


K 


(NuU 

D„, 


(B38) 


The Nusselt and Stanton numbers are calculated from the 
specified or calculated heat transfer coefficient as follows: 


h c D h 

Nu = 

k 

(B33) 

Nu 

St = ■. 

(Re)(Pr) 

(B34) 


As for the pressure drop correlations, CPF provides only 
a limited number (six) of heat transfer correlations for the flow 
path intervals: (1) user-specified heat transfer coefficient, 
(2) trip-strip correlation, (3) pin-fin correlation, (4) finned- 
passage correlation, (5) plain-passage correlation, and (6) input 
Colburn 7- factor curve. There are many other heat transfer 
correlations in the open literature that the user may wish to 
incorporate into CPF (most for a limited and quite specific 
range of geometries). Although the existing heat transfer 
correlations in CPF have a strictly one-dimensional character, 
actual internal heat transfer in turbomachinery shows strong 
three-dimensional effects (such as variations in heat transfer 
between the passage’s leading and trailing surfaces). Such 
variations are accounted for, in a limited sense, by the 
expanded number of correlations used and discussed in 
reference 5. 

For trip strips (roughened tubes, IOPTHT = 1) a correlation 
of reference 6 is used. 


For plain passages the heat transfer coefficient is evaluated 
from an internal correlation (IOPTHT = 4) or from an 
input 7-factor curve (IOPTHT = 5). The internal correlation 
for laminar flow (Re < 2300) is 

h, = 4.226 J- ( Re) 0 °'(Pr)° 333 (B39) 

D h 

This formulation is not of the usual form found in the open 
literature. It was determined from a 7-factor curve fit of 
unpublished engine company data for a rectangular duct with 
a width to height ratio of 3. The internal correlation for 
turbulent flow (Re > 7000) is 

h, = 0.023 (Re)°' 80 (Pr)° 333 (B40) 

D h 

and the internal correlation for transition flow (2300 < Re 
< 7000) is 


K = HL2300 T 


HT7000 - HL2300 


1 


2.951 


1 — exp 


—2.951 (Re — 2300) 
7000 - 2300 



L 

2 


0.28 

4.5(e ,+ ) (Pr) 057 


0.95 



(B35) 


This equation is valid for 10 < p/e < 40 and e + > 35, 
where is given by equation (B26). The heat transfer 


(B41) 


where HL2300 and HT7000 are equations (B39) and (B40) 
evaluated at Re equal to 2300 and 7000, respectively. Note 
that the exact formulation of the laminar and transition 
equations may not be of great importance, since any practical 
cooling design will usually have turbulent flow. 

If the heat transfer coefficient is obtained from an input 
Colburn 7-factor curve (IOPTHT = 5), it is calculated from 
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(B42) 


The Rayleigh and Rossby numbers are given by 


h, = 


JpVC r 

(Pr)"' 6667 

r 


where (Pr)/r is the Prandtl number evaluated at the film 
temperature 


7> = 0.5(7'+ T w ) (B43) 


Ra = 


rN 2 (V)(T K -f)Dl (Pr) 


(f)’ 


(B50) 


For finned passages (IOPTHT - 3) the heat transfer 
coefficient is evaluated from 



(B51) 



m + Tf 


2 tanh (<p f L f ) 
j j + m 

<p f 


where 


Note that equations (B48) and (B49) are valid only for purely 
(B44) radial flow and when the fluid is cooling the wall. If the fluid 
heats the wall (T > 7*,), an appropriate warning message is 
printed out. 


*/ = 


2ft, ,p P 

kTf 



(B45) 


and /? ( pp is the heat transfer coefficient for a plain passage 
obtained from equation (B39). The formulation of equation 
(B44) requires a different heat transfer area to get the proper 
coolant heat pickup (in eqs. (B8) and (B9)). For all heat 
transfer options other than finned passages (IOPTHT ^ 3), 
the effective wetted perimeter is 


Subroutine BYPSFL 

Subroutine BYPSFL (bypass flow) uses simplified versions 
of equations (B8) and (B9) to determine the flow through the 
bypass channel in one marching step. The simplifications are 
the deletion of all heat transfer terms and the use of constant 
physical properties throughout (evaluated at the average static 
temperature of interval NBPI if IOPPRP=l). These 
assumptions require that the bypass channel be kept short and 
that the change in flow area from inlet to outlet be kept 
reasonable. 


P = — (B46) 

For finned passages (IOPTHT = 3) the wetted perimeter is 
based on surface length without the fins, or 

P ~ 2 (Number of fins across the interval)(Fin spacing) 

(B47) 


Subroutine TIPIMP 

Subroutine TIPIMP (tip impingement) calculates the flow 
through an impingement tip cap. The associated flow stations 
are shown in figure 7. Station 2 represents the inlet holes, 
station 3 the impingement chamber at the exit of the inlet holes, 
and station 4 the impingement chamber at the entrance of the 
exit holes. Station 5 represents the exit holes. The chamber 
area is assumed to be much greater than the inlet or exit hole 
areas, giving low throughflow velocity, so that 


Note that this formulation assumes heat soaking into the fluid 
from two sides, as shown in figure 6. 

The correction factors for heat transfer coefficient due to 
rotation (HCCORC(NI)) are from references 9 and 10. For 
radial outflow (ref. 9) 


HCCORC = 


Ra 


(Re)~ 


- 0.186 


(Rs) 


0.333 


(B48) 


and for radial inflow (ref. 10) 


HCCORC =1.72 


Ra 

(Re) 2 


0.112 

(Rs) - 0083 


(B49) 


Pi - Pi (B52) 

From the given inlet and exit conditions and an assumed 
chamber pressure /? 3 , an inflow and an outflow are 
calculated. The chamber pressure is then adjusted in successive 
iterations until the inflow and outflows are equal. Centrifugal 
pumping effects are included from station 3 to station 4 but 
not through the inlet or exit holes. Heat pickup is calculated 
only at the impingement surface (station 4). Heat pickup in 
the inlet and exit holes is neglected. 

Tip-cap inlet flow calculations. — The inlet total conditions 
at station 2 are assumed to be equal to the static conditions 
in the tip-cap interval, as shown by equations (B22) and 
(B23). Inflow static pressure p 2 is assumed to equal the 
chamber total pressure p 2 . This, together with equation (B52), 
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Pi — P 3 — P?> 

For a Fixed chamber pressure p\ an initial guess is made for 
inlet static temperature, and a new temperature is calculated 
from 


2yRg<Ti 


1- * 


* imp = 0.63 — (Re) ; 


x exp — 1 .27 


0 50 /d 2 \ 060 


ti = t 7 + 


050 /D 2 \ 120 


Wir.hp ~ 




(7-l)/7"l 0.50 


7-1 , 

1 +- M\ 

2 


This new temperature is compared with the previously used 
value, and the iteration continues until the calculated 
temperature converges. The total inflow is obtained from 


A 2 total inlet hole area 

D 2 inlet hole diameter 

(SP) 2 inlet hole spacing 

L imp impingement distance 

Z) 4 impingement cavity diameter 

Tip-cap outlet flow calculations .— Outlet flow is calculated 
similarly to the inlet flow. Total exit temperature and pressure 
at station 5 are assumed to be equal to the total conditions at 
station 4, and the exit static pressure is assumed to be equal 
to the specified backpressure. With an initial guess for exit 
static temperature, the calculations proceed as follows: 


w, = CoAy V 2 

“ RgJi 


The discharge coefficient C D is obtained from reference 1 1 
and is shown in figure 12. 

The effects of rotation on pressure from station 3 to station 4 
are calculated by subroutine BYPSFL. The cross-sectional 
flow area from station 3 to station 4 is taken as the product 
of the input impingement cavity diameter (DIC in fig. 7) and 
the calculated length of the tip-cap flow interval (DX(NTCF) ). 
No heat addition or friction is accounted for, since the flow 
path is short and the radial flow Mach number will be low. 
Note, however, that temperature pickup due to impingement 
is calculated as shown in the next section. 

Tip-cap heat transfer calculations .— The temperature 
increase due to impingement in the chamber is calculated by 
the Chupp and Helms equation for jets impinging into a 
concave cavity (ref. 12). 


2t RgXs 


( 7 -D/ 7 I 7 0,50 


yRgjJ- 


( 7 “ I V 7 I 0.50 


The calculated temperature is compared with the previously 
used value, and the calculations are repeated until convergence 
is reached. The exit flow is then calculated from 


(Re) 2 = 


h ? 5 = Qyls F s 

RgJs 



.<1 .6 
MACH NUMBER, M 


Figure 12.— Tip-cap flow discharge coefficient. 


The same discharge coefficient used for the inflow is used for 
the outflow (fig. 12). 

Subroutine TPLOSS 

Subroutine TPLOSS (total pressure loss) calculates the total 
pressure loss due to a sudden expansion or contraction by the 
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resistance coefficients AT/ for pipes (ref. 13). These resistance 
coefficients (or total pressure loss coefficients) are based on 
the velocity in the smaller pipe and are defined by 




Ap' 



smaller pipe 


(B65) 


For sudden expansion the total pressure loss coefficient can 
be expressed analytically as 


A \\ 2 

1 L (B66) 

A 2/ 

where A\ is the inlet area and A 2 is the exit area. The curve 
for sudden contraction was curve fit by the following 
expressions: 

/ A \°50 

*=(£) <b671 



K T = -0.41468 X 3 + 0.07014 X 2 - 0.25946 X + 0.50 

(B68) 


for 0.0 < X< 0.80, and 


K r = 3.64583 X 2 - 7.1875 X + 3.54167 (B69) 


Equation (B71) is a transcendental equation and is solved by 
Newton's method. 

Roughened passage.— A passage js considered to be rough 
if the specified relative roughness Die (input parameter DOE) 
is less than l.Ox 10 6 . For this case reference 7 suggests the 
following formulas for friction factor: 


/ = — for Re < 2300 (laminar flow) (B73) 
Re 

For fully developed turbulent flow the friction factor ex- 
pression is not a function of Reynolds number and is given by 

1 (D\ 

~jr = 1 .737 In - + 2.28 (B74) 

v / w 

This formula is valid for 


< 0.01 (B75) 

For the transition region the friction factor is a function of 
both Reynolds number and relative roughness and is given by 



for 0.80 < X < 1 .0: where A\ is again the inlet area and A 2 , 
the exit area. 


1 / D\ 

— — - 1.737 In - ) + 2.28 
V/ W 




Subroutine WWWFF 

Subroutine WWWFF (Welty, Wicks, and Wilson friction 
factors) establishes friction factors for both smooth pipes and 
rough pipes as shown in reference 7. 

Smooth passage .— A passage is considered to be smooth 
if the specified relative roughness Die (input parameter DOE) 
is greater than 1 .Ox 10 6 . Reference 7 suggests the following 
formulas for friction factor: 


f — — for Re < 2300 (laminar flow) (B70) 
Re 


- 1.737 In 


4.67 


(Re)v7 


The equation is valid for 


(B76) 



(Re)Vy* 


> 0.01 


(B77) 


and is again solved by Newton’s method. 


-L = 1.737 ln[ (Re)V/| - 0.40 

,i f 

for Re > 3000 (fully turbulent flow) (B7 1 ) 

For the Reynolds number range 2300 to 3000 the following 
equation is assumed: 

„ ln(Re) -7.2710 


Subroutine AIRPRP 

Subroutine AIRPRP (air properties) calculates physical 
properties for unvitiated air from ambient to stoichiometric 
temperatures. The physical properties are specific heat ratio, 
specific heat, viscosity, and thermal conductivity. They are 
evaluated from curve fits of the charts of reference 14. Below 
1650 °C (3000 °F) the calculated values are independent of 
pressure; above that temperature the generated values are valid 
only for a pressure of 3 atmospheres. 


67.5144 



Subroutine SPLFIT 


Subroutine SPLFIT (spline fit) generates a spline curve fit 
from each input set of tabular data. The tabular data set must 
consist of at least 3 but not more than 25 points. The curve- 
fitting procedure requires the slopes at the end points. These 
slopes are calculated from the first two and last two data points. 


For this reason these points must be chosen such that a straight 
line between them gives a good approximation to the slope 
of the curve at the end points. If the calling program calls for 
a value at an x location outside the range of the input table, 
the value at the nearest end point is used. 
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Appendix C 
Input Variables 


Input 

Description and units 

Variable 

Default 

parameter 1 


type b 

value 

Overall parameters 

TITLE 




I ECHO 

Option for input echo: 


I 

1 


0— no input echo 





1— input echo 




ISLOUT 

Option for printing output at slice locations: 



0 


0— print only at original input nodes 





1— print at all slices 




IUNITS 

Units for input variables: 



0 


0— English units 





1— SI units 




IOPPRP 

Option for physical properties specification: 



1 


0— user specifies G, CP, XMU, and XK 




1— properties obtained from subroutine 




AIRPRP 



G 

Specific heat ratio 0 


* 


CP 

Specific heat at constant pressure, 1 





cal/(g K); Btu/(lbm °RJ 




XMU 

Viscosity, 0 g/(cm s); lbm/(ft s) 




XK 

Thermal conductivity/ cal/(cm s K); 





Btu/(ft hr °R) 




RG 

Gas constant, J/(kg K); (ft lbf)/(lbm °R) 




RPM 

Rotative speed, rpm 

i 



NI 

Number of intervals in flow passage 

i 


ZIN 

Inlet z location, cm; in. 

R 


RIN 

Inlet r location, cm; in. 




YIN 

Inlet y location, cm; in. 




AIN 

Inlet free flow area, cm 2 ; in. 2 




DHIN 

Inlet hydraulic diameter, cm; in. 




IOPTW 

Option for coolant flow rate: 

I 



0— user specifies coolant flow rate (WIN) 




1— program calculates coolant flow rate 




from specified inlet and exit pressures 




(PTSP and POUT) 



WIN 

Coolant inlet flow r rate, d kg/s; Ibm/s 

R 


PTSP 

Inlet total supply pressure, N/cm 2 ; psa 




TTSP 

Inlet total supply temperature, K; °F 




KTIN 

Inlet total pressure loss coefficient 



0.50 

POUT 

Outlet pressure, 0 N/cm 2 ; psia 

V 



Input curves 

NPFF2 

Number of points describing Fanning friction 

I 



factor curve for pin fins (3 < NPFF2 < 25) 



IXFF2 

X axis coordinate indicator: 

I 

0 


0— Cartesian 




1— logarithmic (base 10) 



RE2 

Values of Reynolds number at NPFF2 input 

R(NPFF2) 



points in ascending order 



IYFF2 

v axis coordinate indicator: 

I 

0 


0— Cartesian 




1— logarithmic (base 10) 



FFF2 

Fanning friction factor values at input points 

R(NPFF2) 


NPFF3 

Number of points describing Fanning friction 

I 



factor curve for finned passage 



. 

(3 < NPFF3 < 25) 






Input 

Description and units 

Variable 

Default 

parameter 3 


type 1. 

value 


Input curves (concluded) 


IXFF3 

x axis coordinate indicator: 

I 

0 


0— Cartesian 




1— logarithmic (base 10) 

R(NPFF3) 


RE3 

Values of Reynolds number at NPFF3 input 



points in ascending order 

I 

0 

IYFF3 

v axis coordinate indicator: 


0— Cartesian 




l — logarithmic (base 10) 

R(NPFF3) 


FFF3 

Fanning friction factor values at input points 


NPFF4 

Number of points describing Fanning friction 

I 



factor curve for plain passage 
(3 < NPFF4 < 25) 


0 

1XFF4 

x axis coordinate indicator: 
0 — Cartesian 

1 



1— logarithmic (base 10) 

R(NPFF4) 


RE4 

Values of Reynolds number at NPFF4 input 



points in ascending order 

I 

0 

IYFF4 

y axis coordinate indicator: 


0— Cartesian 




1 — logarithmic (base 10) 

R(NPFF4) 


FFF4 

Fanning friction factor values at input points 


NPCJF 

Number of points describing Colburn 7-factor 

I 

— 


curve (3 < NPCJF s 25), 
J = MPrr’/pFC,, 

I 

0 

IXCJF 

x axis coordinate indicator: 


0— Cartesian 




1— logarithmic (base 10) 



RECJF 

Values of Reynolds numbers at NPCJF input 

R(NPCJF) 



points in ascending order 


0 

IYCJF 

v axis coordinate indicator: 

I 


0— Cartesian 




1— logarithmic 

R(NPCJF) 


CJF 

Colburn 7- factor values at input points 


Nodal and interval parameters 


Z nd 

z location, cm; in. 

R(NI) 


R",/ 

r location, cm; in. 

R(NI) 


ytui 

v location, cm; in. 

R(NI) 


NS LI' 

Number of slices in each interval 

KNI) 

1 

IOPTFP 

Option for fake interval specification: 

KNI) 

0 


0— regular interval 




1— fake interval 



IOPTHT' 

Option for interval heat transfer calculations: 

I(NI) 

0 


0— user specifies h via input parameter HC; 




for this option the user must also specify 
friction factor via input parameter F 




1 — trip strips (turbulator correlation of 
ref. 6); user must input EDD and PDE 

2— pin fins (Van Fossen correlation. 

ref. 8); user must input AV and DHV 

3— finned passage; user must input NFNP, 




THFN, SPFN, HHTFN, and KMET 




4— plain passage (internal correlation) 

5— user specifies Colburn 7-factor curve 




via input parameters NPCJF, IXCJF, 
RECJF. IYCJF, and CJF 





input 

parameter 1 

Description and units 

Variable 

type 11 

Default 

value 

Nodal and interval parameters (concluded) 

ioptff 

FCOR' 

HCCOR' 

HTACOR' 

IOPTCR' 

EDIT 

PDF.' 

DOB' 

KTSGMT' 

NFNP' 
THEN' 
SPFN' 
HHTFN' 
KM BP 

A" <y 

AV' 

HC' 

DH' i(/ 

DHV' 

F' 

TW' 

Q 1 

Option for Fanning friction factor 
specification: 

0— user specifies / for each interval via 
input parameter F 

1— trip strips (turbulator correlation of 
ref. 6); friction factors generated from 
input parameters EDD and PDE 

2— obtained from input curves specified 
by input parameters NPFF2. etc., for 
pin fins; input parameters NPFF3, etc., 
for finned passage; input parameters 
NPFF4, etc., for plain passage 

3 — obtained from reference 7 correlation 
for roughened plain passage (must 
specify DOE) 

4 — obtained from specified total pressure 
loss coefficient KTSGMT 

Friction factor correction 
Heat transfer coefficient correction 
Heat transfer area correction 
Option for calculating rotational heat transfer 
coefficient correction: 

0 no correction 

1— correction (specify only for RPM > 0.0) 
Relative roughness of trip strips’ 

Relative rib spacing of trip strips’ 
Diameter-to-roughness ratio* 1 in reference 7 
(roughened-passage friction factor 
correlation) 

Total pressure loss coefficient (may need to 
correct for path length via HTACOR) 
Number of finned passages in interval’ 1 
Fin thickness, 11 cm; in. 
bin center-to-center spacing. 1 ' cm; in. 
Half-height of fin. 11 cm; in. 

Metal thermal conductivity,’ 1 cal/(cm s K); 

Btu (ft hr °R) 

Cross-sectional area based on minimum free 
flow area, cm 2 ; in. 2 

Cross-sectional area based on open volume, 1 
V/L, enr; in. 2 

Coolant heat transfer coefficient, 1 
cal/(cnr s K); Btu/ft : hr °R) 

Hydraulic diameter based on minimum free 
flow area (D h = 44/P), cm; in. 

Hydraulic diameter based on open volume, 1 
4 V/ 5, cm; in. 

Fanning friction factor 11 

Inside wall temperature, K; °F 

Heat generation, J/kg; Btu/lbm (see appendix E) 

KNI) 

R(NI) 
R(NI) 
R(N1 ) 
KNI) 

R(NI) 

KNI) 

R ( N I ) 

T 

0 

1.0 

1.0 

1.0 

0 

0.0 

().() 

0.0 

Parameters specifying locations of tip-cap and bypas; 

i now 


NTCF 

NBPI 

NBPO 

Interval from which tip-cap flow is specified 
Interval containing start of bypass channel 
Interval containing exit of bypass channel 

i 

i 

i 

— 
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Input 

parameter* 

Description and units 

Variable 

type b 

Default 

value 


Tip-cap parameters for interval NTCF 


lOPHTT 

Option for tip-cap How heat transfer 

1 

0 


calculations: 




0 — no heat transfer calculations (How only) 




1— heat transfer based on reference 12 




(correlation for impinging jets into 




concave cavity) 



NHI 

Number of tip-cap inlet holes 

I 

— 

DHI 

Diameter of tip-cap inlet holes, cm; in. 

R 


NHO 

Number of tip-cap outlet holes 

1 



DHO 

Diameter of tip-cap outlet holes, cm; in. 

R 


PBP 

Tip-cap outlet hole backpressure, N/enr: 





psia 




RIH 

Radius at inlet holes (must be greater than 





R(NTCF)), cm; in. 




R1PL 

Radius at impingement plane (must be greater 





than RIH), cm; in. 




SPIH 

Spacing of impingement holes, cm; in. 




ALIMP 

Angle between main-stream flow and 



90.0 


impingement flow. 1 deg 




DIG 

Diameter of impingement cavity, cm; in. 




AIMP 

Impingement cavity heat transfer area, 





cm 2 ; in.~ 




TWIMP 

Impingement cavity wall temperature, K; F 





Bypass channel parameters 



DLBP 

Channel length, cm; in. 

R 

— 

DHBP 

Channel hydraulic diameter, cm; in. 




FFBP 

Channel Fanning friction factor 




AINBP 

Inlet area, cm 2 ; in." 




PBPRI 

Pressure recovery at inlet 



0.0 


(0.0 < PBPRI < 1.0) 




RINBP 

Inlet radius, cm; in. 



— 

ALBPI 

Angle that describes component of exiting 



90.0 


flow in main-stream direction in interval 





NBPI, 1 deg 




AOUTBP 

Outlet area, cm 2 ; in." 




ROUTBP 

Outlet radius, cm; in. 




ALBPO 

Angle that describes component of entering 



90.0 


flow in main-stream direction in interval 





NBPO (tig. 3), deg 





a Supcrscnpl i denotes an interval value; supcrseripi 
tul denotes a nodal value. 

h I denotes unsuhseripted inleper value; l< I denotes 
subscripted integer value. 

R denotes unsuhseripted real value; R( ) denotes 
subscripted real value. 
c lnput only if IOPPRP - 0 
‘'input only it IOPTW - 0 


c lnpm only il IOPTW - I 
'input only if IOPTHT - I 
^ Input only if IOPTFK J 
'‘input only if IOPTHT - T 
'input only if 1POTHT - 2 
'input only if lOPTH I' - 0 
14 Input only if IOPTI I - t) 

'Dummy variable (not presently used! 


ORIGINAL PAGE IS 
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Appendix D 

Derivation of One-Dimensional Momentum and Energy Equations 


Equations (B8) and (B9) are derived in this appendix. They 
represent formulations of the one-dimensional momentum and 
energy equations, respectively, which can be integrated along 
the flow path. 


The centrifugal pumping force component acting on the control 
volume and directed along the flow-path axis is expressed as 
mass times acceleration by 


Momentum Equation 

Consider the flow of fluid through a control volume in a 
flow path rotating with speed N about an axis as shown in the 
sketch. The distance along the flow path is a, and the distance 
from the rotation axis is r. 



Fluid is injected at a rate dw at velocity F inj , having a velocity 
component in the main-stream flow direction of V injx . The 
momentum equation states 


Change in momentum = ^ Forces (pressure, friction, 
flow rate (out-in) and body) 


The friction force is expressed by the Fanning friction factor as 


P 



(D4) 


where the angle a is as shown in the preceding sketch. Since 


cos a = — (D5) 

dx 

and after neglecting second-order terms (products of differ- 
entials), the centrifugal force component becomes 

rN 2 pA dr (D6) 

The forces acting on the control volume are represented in 
the following sketch: 



The body force due to the Coriolis acceleration will always 
act normal to the flow path and makes no contribution along 
the flow direction. The momentum equation thus becomes 


- pV 2 fP dx 


From the definition of hydraulic diameter 


4A 

Du — — 


the friction force becomes 


- pV 2 f— dx 
2 D h 


(w + dw)(V + dV) -\wV+ V miX dw 


(Dl) 


— pA + (p + ^IpjdA — (p + dp) (A + dA) 

pV 2 f4A 


2D„ 


dx + rN'pA dr 


(D7) 


(D2) 


Expanding and dropping second-order terms results in 


wdV + Vdw - V in : x dw = -A dp - 


P V 2 f4A 


dx 


(D3) 


2 D h 

+ rN 2 pA dr (D8) 


28 



Dividing by dx and combining wdV and Vdw gives 


The momentum equation can then be written as 


d . dw dp p V 2 fAA , J dr ( 2wRg c T \ dw ^ w 2 Rg, dT w^RgJ dA 

— (wK) - J'V — - ~ A ^ 2 D h +r P dx \ pA m,x ) dx pA dx pA 2 dx 


dx 

(D9) 

The left-most term is expanded by using the perfect-gas 
law and 


Thus, 


wV = 


pA 


J 

pA 


and 


A (wv) 

dx 


2wRg c T dw . , 
pA dx 

l dT T dA T dp 

<A dx pA 2 dx Ap 2 dx 


(DIO) 


(Dll) 


(D12) 


w2r xJ dP = _ A ±__ f> y2 f 4A + rN 2 pA ± 

Ap 2 dx dx 2 D h dx 

(D13) 


Solving for dp/dx gives 


dp f A _ w 2 Rg, T 
dx \ Ap 1 


-2wRf!,T \dw 

pA inj 7 dx 

w 2 Rg r \ dT / w 2 Rg c T \ dA 
pA ) dx \ pA 2 ) dx 


p V 2 f4A dr 

f* rN pA — 

2 D h dx 


(D 1 4) 


Using the perfect-gas law and equation (DIO) in the friction 
and pumping terms and dividing through by A gives the tinal 
form of the one-dimensional momentum equation. 


dp 

dx 




inj,x 


pA 


(i) R ± dT | (jW T_ M , Mp dr 

dx p dx pA dx 2 pD h RgX dx 




(D15) 


Energy Equation 

Consider the flow of fluid through the control volume shown 
in the following sketch: 


states that 


Change in energy = ^ Heat fluxes by external sources 
flow rate (out-in) 


d w 



w + dw 


v + dV, h + dh 


Fluid is injected at a rate dw w ith a velocity V my The energy 
in a gas stream is given by h + K 2 /2 and the energy equation 


Assuming that the added heat is due to convection (via the 
convective heat transfer coefficient h c ), the energy equation 
can be expressed as 


(w + dw) 


h + dh + 


(V+dVy 


v l \ 

h+ — | vv 


Vjn 

2 


h ini + ^ )dw = hf dx(T u . - r aw ) 


(D16) 
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After expanding the left-hand side and neglecting second-order 
terms, this reduces to 


wC — + wRg > T ( wR z< :r 

'' dx pA dx\ pA 


vv (III 4* wV (IV 4- (h — h- m i)dw 4- — 


V 2 vl 


dw 


,nj ' '2 2 

= h,P dx(T w - T a J (D17) 


Dividing by dx and using 


h = C p T 


wRgJ 


Ap 


the energy equation can be arranged to give 


1 / wRg c T 

r T - C T A — I — 

p pmy mi 2\ pA 


- vi 


mj 


dw 

dx 


= h c P(T w — r aw ) 



The second term can be expanded to give 

(018) 

d . 

( wRg c T\ 

Rg c T dw wRg c dT 
— 4* 


dx 

l P A ) 

pA dx pA dx 

(D19) 



wRg c T dA wRg c T dp 


(D20) 


— (D21) 


pA dx p l A dx 

Substituting this back into the energy equation and collecting 
the terms containing dT/dx gives 


— ( w c + W ~ R & T > v /? g - 

dx V p pA pA 


RgJ CT+C l (wRg c TY 1 2 

pA pA p /M " J ,nj 2 \ pA ) 2 ‘ nj 


dw 

dx 


+ 


r w 2 Rg c T wRg,.T 

dA 

_L 

£ 

dp_ 

1 

tsj 

1 

dx 

K) 

1 

dx 


r.w) 


(D22) 


Dividing through by wC p and using the definition of hydraulic 
diameter (eq. (D2)) gives the final form of the one-dimensional 
energy equation. 


dT _ 1 
dx~ CT 


r 

L ^ r 


C T 

'“■/Mnj J in j 

C p w 


pA 

vi 


w fRg c T\ + T + w ( R gl T 


w 2 C„ \ pA 


mj 


2hC 


dw 

dx 


+ 


+ 


J ( wR *x 

CpA \ pA 

4h,A{T w ~ T m ) 


dA 

— + 

1 

( wRgJ\ 2 

dx 


\ PA ) . 


wC p D h 


where 


CT= 1 4 


Marching Equations 




dp 

dx 


Equations (D15) and (D23) are two simultaneous equations 
in the variables dp/dx and dT/dx . In order to march the 


solution along the flow path, they must be solved for the 
individual variables. After substituting the energy equation into 
the momentum equation (replacing dT/dx) and performing 
considerable algebraic manipulation, the momentum equation 
takes the form shown as equation (B8). 


dp _ 1 
dx ” C2 


I 

r 

1 

) C1 

i 




(D23) 

1 

+ < 

< 



'N 

wRgjV i 


1 pA 

CpA \ 

pA ) \ 

J 

(D24) 


. 2 




4JM5) 

N 2 rp dr 

4h,AC\ 



2pD h 

RgJ d x 

w 


dA 

dx 


wC,,D h 


(D25) 
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where 


Cl = - 

P 



(D26) 


C2 





(D27) 


This value for dp/dx is then substituted into the energy 
equation (D23) to evaluate dT/dx, and the solution is marched 
from node to node as explained in appendix B. 
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Appendix E 

Feature of Computer Code for Predicting Coolant Flow and 
Heat Transfer With Specified Heat Generation 

Teresa R. Kline 


This appendix describes a feature of the computer code CPF 
that allows the direct specification of heat generation in a given 
flow interval. With this feature Rayleigh line flow problems 
can be solved directly. The appendix provides derivations of 
the equations, shows input and output, and demonstrates the 
use and accuracy ol the feature via an example problem. 

The code originally had no provision for direct heat generation 
(specifying a desired quantity of energy addition or subtraction 
in a given flow interval). Heat addition or subtraction was 
accounted for through a specified or calculated heat transfer 
coefficient, a heat transfer area, and the temperature difference 
between the wall and the fluid. This formulation made the 
quantity of heat addition or subtraction dependent upon the 
flow path length. The main text demonstrates the accuracy of 
CPF via example problems for adiabatic flow, Fanno line flow, 
and Rayleigh line flow. Since Rayleigh line flow involves the 
direct addition or subtraction of specified quantities of heat 
(independent of path length), the Rayleigh line flow example 
problem could not be solved directly. Instead, fixed values 
of wall temperature and heat transfer coefficient were assumed, 
and the path length (heat transfer area) was varied until the 
desired exit Mach number was reached. The calculated exit 
temperature and pressure were then compared with the closed- 
form Rayleigh line solution values. 

This appendix describes a feature of CPF in which a direct 
heat generation term (energy per unit mass) added to the energy 
equation allows the direct solution of the Rayleigh line flow 
problem. This added feature can also be used to account for 
heat of vaporization (CPF allows the user to specify gas fluid 
properties other than those for air). In addition, the revised 
code can be used to calculate and compare results from heat 
transfer tests, where known quantities of heat are often added 
(e.g., via electric heater strips). 

This appendix also describes the input and output of the 
modified code. An illustrative example is given to show the 
effects of both positive and negative heat generation. The 
formulation of the energy and marching equations is given with 
the heat generation term added. CPF is available through 
the Computer Software Management Information Center 
(COSMIC), University of Georgia, Athens, Georgia, 30602. 

Input 

The modified code contains one new input variable, the heat 
generation per unit mass, Q(I). It is an interval variable and 
must be specified for each interval. The units are British 
thermal units per pound mass (Btu/lbm) in the U.S. customary 
system or joules per kilogram (J/kg) in the International 
System of Units (SI). The heating term is defaulted as 0.0. 


If Q = 0.0, we are simply left with the original code. The input 
value for heat generation is echoed in the output listing, where 
it is located immediately before the geometric output. Correct 
units (U.S. or SI) are included. 

Derivation of Energy Equation 

The momentum equation is the same as that given in the 
main text, but the energy equation has an added heat generation 
term. The energy equation is derived as follows: 


dw 



w + dw 


V + dV, h + dh 


As shown in the sketch, the change in energy flow rate 
equals the sum of the heat fluxes by external sources: 


(w + dw) | ( h + dh) + 


(V + dvy 


w ( h + 


V 2 


dw^h mj 4- ^ J ^ — h c P (T w — F aw ) dx + wq 


Expanding and neglecting second-order terms gives 


w dh + wVdV+h dw + } L*!L _ h dw _ Yk dw 

2 ,nj 2 


— h c P ( T w — T aw ) dx T wq 


Grouping terms gives 

vv dh + wV dV + (h - h m] )dw + 




dw 


= h c P (T w - T aw ) dx + wq 
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Substituting h = C p T and V = wRg c T/pA gives 


Solving equation (El) for dT/dx yields 


w'-Rg.T ( wRg c T 

wC p dT + — d[ — 

P pA \ P A 


+ ( 1 CpT C p ;\r\)T\r\) + 2 p 2 A 2 




= h,P{ T w - 7 aw ) dx + wq 


. w 3 R 2 g*T\ dT_f 

VVC., H 1 — l H- 


/r/4 2 ) dx 


w 2 /? 2 g^ , H |2 /? 2 g?:r 2 \ ^ 

2/>vl 2 2 p-A 2 ) dx 

w 2 R 2 g;T 2 dp w'R 2 g;T 2 dA 
p } A 2 dx p 2 A 3 dx 


Dividing by dx yields 


+ h r p(r K - t. m ) + 


. ^ dT w 2 Rg c T d fwRg e T 

wC n 1 I 

P dx pA dx \ pA 


+ f C d T - C„i niTini + 


w 2 R 2 g 2 .T 2 Vi.\ dw 


2 p^A* 2 } dx 


Dividing by wC p and substituting P — 4 A/D h gives 


w 2 R 2 g 2 .T \ dT _ / T + CpjnjT'inj _ wR~gfT~ 
P 2 A 2 C p ) dx ~ \ w CpW 2 C pP 2 A 2 


~ hcP\Jw T w ) + ^ 


Expanding the second term gives 


d_ ( wRg c T 
dx\pA 


Rg c T dw wRg c dT 

pA dx pA dx 


wRg c T dp wRg c T dA 

p 2 A dx pA 2 dx 


VU wR 2 g 2 T 2 \ dw 
2 C p w C p p 2 A 2 ) dx 

w 2 R 2 g;J 2 cip w 2 R 2 glT 2 dA 

C p p } A 2 dx C p p~A* dx 


h<AA(T w - r aw) q 
1 C p wD h C p dx 


From the momentum equation 


2 W**' T , ^ 


S Rg ' dT | 


j) R *x dA wx-nt 

pA dx 2pD h 


N 2 rp dr 
+ RgJ dx 


i) s *' r 


From the energy equation 


_ T C,,, inj T, nj _ (Rg£\ 2 + _ A (**£)'] + _L (***£)' — — ( wRgJ \ — + /if4M(r ': ~ ^ 

»• + Cm 2 C„ \pA ) 2 Cpv C r \ pA ) dx C„p \ pA ) dx C,/t \ pA J dx C p wD k 


O Rl & 
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Substituting dT/dx into dp/dx gives 



1C \ pA 


3VV (Rgj\ | T Cp jnj 7] n j 

w C fJ w 


V 2 • 

r inj 

2wC, 


7>J 


mj 




pA 

■> 

vM5)’ 


Cl / wRg c T 
C p A \ pA 


dw 

dx 


dA 

dx 


N 2 rp dr 
Rg c Tdx 


pD h 

4M(7’ H . - r aw )Cl C\q 


w CpPh 


Cp dx 


where 


(station 1 to 2 ) and the cooled interval (station 2 to 3 ) shown 
in the following sketch: 



INTERVAL 2-3 IS A MIRROR IMAGE OF INTERVAL 1-2 


The problem can be stated as follows: 

*7in out 

The intervals are of equal length and diameter. Begin with 
set initial conditions: P 0 j, r 01 , A, 7 , D, vv, and A/,. Add heat, 
< 7 in = 1 329 288 J/kg, between station 1 and station 2. Take 
away heat, q out - - q m , between station 2 and station 3 to see 
if this gives (within computational error) the original input 
values. 

From Shapiro (ref. 15) the closed-form solution is as 
follows: 

Solve for T 02 : 


Cl = 


© 


1) **- 


5) R2g ' J 

P 2 Cp 


C 2 = 1 




Cl wRgJ 


pC \ pA 


^in ~ C p ( T()2 T { ) |) 


Iterate for M 2 : 

r 02 _ m \ (1 +7^t) 2 1 

7 !ii 0 + 7 A/T) 2 7-1 , 

V V 1 + 

2 


Solve for /? 2 , /? U2 , and T 2 : 


This value of dp/dx is then substituted into the energy equation 


P2 _ 

1 + yM] 

to give dT/dx. The solution is then marched from node to node 


P\ 

i +yM\ 

as described in the main text. 




Example Problem 

Poi_P 2 ( 

1 

2 


Po I P\ / 

, 7 - 

1 _ _•> \ *y /CV ~ 1 1 

The example given here is a reworked and expanded 

( 

1 + - — 

-Mt) 

example 3 from the main text. It is the Rayleigh line flow 


z 


problem, which could not be solved directly with the original 




code. It involves adding heat (1329.288 kJ/kg) between stations 

t 2 

m\ ( 

1 + 7 Mt) : 

1 and 2 and subtracting the same amount of heat between 

T 

■ jf2 / 

/ 

stations 2 and 3. In the main text only heat addition is shown. 

T 1 


] + 7 M \\ 2 


This example doubles the length of the duct, demonstrates that 
heat subtraction is also allowed by the code, and shows how 
closely the end solution matches the input conditions when an 
equal amount of heat is added and subtracted. For this example 
equal numbers of slices were used for the heated interval 


Results are shown in table E-I for 50, 100, and 125 slices 
in each interval. Improvement is evident between 50 and 100 
slices, but between 100 and 125 slices improvement is so slight 
that additional slices cannot be justified. 
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TABLE E-L— COMPARISON OF RESULTS 


(a) Heating in a duct 


Conditions at 
duct midpoint 

Closed form 

50 Slices 

100 Slices 

125 Slices 

7V 

1625 

1625 

1625 

1625 

T 2 

1461 

1461 

1461 

1461 

Pi 

83.44 

83.44 

83.44 

83.44 

Pi 

57.45 

57.48 

57.47 

57.47 

m 2 

0.750 

0.7496 

0.7497 

0.7497 


(b) Cooling in a duct 


Conditions in 
duct exit 


Closed form 

50 Slices 

100 Slices 

1 

125 Slices 

300 

300 

299.9 

299.9 

297.6 

297.6 

297.5 

297.5 

100 

99.99 

99.99 

99.99 

97.25 

97.24 

97.24 

97.24 

0.2000 

0.2000 

0.2000 

0.2000 


h' 

h 

Pi 

Pi 

Mi 


The input and output for the example problem are as follows: 


EXAMPLE PROBLEM 2 HEAT GENERATION - SI UNITS 
SDATT IECHO=0 , IUNITS=1, IOPPRP=0, RG=287.0, RPM=Q.O, 

G=1.40, CP=0. 23978, XMU=l.BE-4, XK=0.6E-4, 

NI=2, ZIN=0 . 0 , RIN=0 , 0 , YIN*0.0, AIN=10., DHIN=3 . 56825 , 

I 0PTW=0 , WIN=. 787355, 

PTSP = 1 0 0 . , TTSP-300., KTIN=0.0, 

2=10 0, 20 0, R=2*0 . 0 , Y=2#0.0, NSLI=2*50, A=2*10., MC=2«0.0, 
DH=2*3. 56825, F=2*0.00, TW=2*533., Q=1329287 . 692, -1329287.692 SEND 
SDATT NSLI=2*100 &EHD 
SDATT NSLI=2*1 25 SEND 


EXAMPLE PROBLEM 2 HEAT GENERATION - SI UNITS 


SI UNITS 

USER SPECIFIES PHYSICAL PROPERTIES CG, CP, MU AND K) 

COOLANT FLOW RATE IS SPECIFIED 

ROTATIONAL SPEED IS O.OOOE+OO RPM 

OUTPUT IN ORIGINAL NDDE/INTERVAL FRAME OF REFERENCE 


INPUT 

SLICES 

SLICE 

SPECIFIED HEAT 

NODE 

IN 

NO. 

IN INTERVAL 

NO . 

INTERVAL 


(J/KG) 

1 

50 

50 

0 . 13293E+07 

2 

50 

100 

-0. 13293E+07 


I 

OPTIONS 

( IOPT) 

z 

R 

Y 

A 

DH 

TW 

DX 

XT OT 

KT 

NSL I 


FI 

HT 

FF 

CR 

(CM) 

(CM) 

(CM) 

<CM**Z) 

(CM) 

<K) 

(CM) 

(CM) 








0 . OOOOE+OO 

0 . OOOOE+OO 

0 .OOOOE+OO 

0 .1000E+02 

0 .3568E+01 




0 . OOOOE+OO 


1 

o 

0 

0 

0 

0 . 1000E+02 

0 . OOOOE+OO 

0 .OOOOE+OO 

0 . 1000E+02 

0 . 3568E+0 1 

0 . 5330E+03 

0 . 1000E+02 

0 . 1000E + 02 

0 . 0000E+ 00 

50 

2 

0 

0 

0 

0 

0 .2Q00E+Q2 

0 .OOOOE+OO 

0 .OOOOE+OO 

0 . 1 Q0QE+02 

0 . 3568E+01 

0 . 5330E+03 

0 . 1000E* 02 

0 . 2000E+02 

0 . OOOOE+OO 

50 


I 

T 

TT 

p 

PT 

W 

V 

M 

GAMMA 

CP 

vise. 


TH. CNDT . 


(K> 

(K) 

(N/CM**2) 

(N/CM#*2) 

(KG/S) 

(M/S) 


CAL/(GRK) GR/( CM- S ) 

CAL/CCM-S-K) 


0 . 2976E+03 

0 . 3000E+03 

0 . 9725E + 02 

0 . I 000E+Q3 

0 . 787AE+ 00 

0 . 6915E+02 

0 .2000E+00 

0 . 1A0E+01 

0 .2A0E+00 

0 . 180E- 

03 

0 . 600E-09 

1 

0 . 1 A61E + 0A 

0 . 1625E+0A 

0 . 57 98E + 02 

0 .8344E+02 

0 .7B7AE+00 

0 . 5793E+03 

0.7496E+00 

0.14DE+01 

0 . 240E+0Q 

0 . 1O0E- 

03 

0 . 600E-04 

2 

0 . 2976E+03 

0 , 3000E+03 

0 .9724E+02 

0 .9999E+02 

0 .787AE+00 

0 . 6915E+02 

0 .2000E+00 

0 . 140E+01 

D . 240E+00 

0 . 180E- 

03 

0 . 600E-0A 


RE 

PR 

ST 

NU 

F 

FCOR 

HC 

HCCOR 

HCCORC 

HTACOR 


TAW 







<CAL/CM**2-$-K) 





( K ) 

1 

0 . 1 561E+07 

0 .7194E+00 

0 .OOOOE+OO 

0 .OOOOE+OO 

0. OOOOE+OO 

0 . 1000E+D1 

0. OOOOE+OO 

1 .0000 

1 . 0000 

l . onoo 

0. 

1595E+04 

2 

0 . 1561E+07 

Q.7194E+00 

0 .OOOOE+OO 

0 .OOOOE+OO 

0. OOOOE+OO 

0 . 1 OOOE+O 1 

0 .OOOOE+OO 

1 .0000 

1 .0000 

1.0000 

0 . 

, 3129E + 03 
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EXAMPLE PROBLEM 2 HEAT GENERATION - SI UNITS 


SI UNITS 

USER SPECIFIES PHYSICAL PROPERTIES (G, CP, MU AND K) 

COOLANT FLOW RATE IS SPECIFIED 

ROTATIONAL SPEED IS O.OOOE+OO RPM 

OUTPUT IN ORIGINAL NODE/INTERVAL FRAME OF REFERENCE 


INPUT SLICES 

SLICE SPECIFIED HEAT 



NODE IN 

NO. IN 

INTERVAL 



NO 

INTERVAL 


( J/ KG) 



1 

100 

100 0. 

13293E+07 



2 

100 

200 -0, 

13293E+07 



I 

OPTIONS ( IOPT ) Z 

R 

Y 



FI HT FF 

CR (CM) (CM) 

(CM) 




0. OOOOE+OO 0. OOOOE + OO O.OOOOE + OI 

1 

0 0 0 

0 0. 1000E+02 0. OOOOE+OO O.OOOOE+OI 

2 

0 0 0 

0 0 . 2000E + 02 0. OOOOE + OO O.OOOOE+OI 

I 

T 

TT 

P 

PT 



(K) 

(K) 

(N/CMXX2) 

(N/CM**2) 



0 . 2976E+03 

0 . 3000E + 03 

0 . 9725E+02 

0 . 1 D00E + 03 

0 

1 

0.1A61E+04 

0 . 1625E + 09 

0 .5747E+02 

0 .83AAE + 02 

0 

2 

0 . 2975E + 03 

0 .2999E + 03 

0.9724E+02 

0 .9999E + 02 

0 


RE 

PR 

ST 

NU 


1 

0. 1561E+07 

0 . 719AE + 00 

0 .OOOOE+OO 

0. OOOOE + OO 

0 

2 

0 . 1561E + 07 

0 . 719AE + 00 

0 . OOOOE+OO 

0 .OOOOE + OO 

0 


A 

CCMX#2) 

0 . 1 000E+02 
0. 1000E + 02 
0 . 1000E+02 


DH 

(CM) 

0 .3568E+01 
0 .3560E+01 
0 .3568E+01 


TW 

(K) 


0 .5330E + Q3 
0 . 5330E+03 


DX 

(CM) 


0. 1000E + 02 
0 . 1Q00E + 02 


XTOT 

(CM) 


0 . 1 0O0E+02 
0 .20Q0E + 02 


0 .OOOOE+OO 
0. OOOOE + OO 100 
0. 0000E+00 100 


w 

(KG/S) 


7074E+00 


V 

(M/S) 

0,691 5E+02 
0 .57A3E + 03 
0 .691AE+02 


GAMMA CP 

CAt /(GR-K) 


vise. 

GR/( CM- S) 


TH. CNDT . 
CAL/( CM-S-K ) 


0 .2000E + 00 
0 . 7497E+00 
0 .2000E + 00 


0. 140E+01 
0 . 1A0E+01 
0 . UiOE+Ql 


HC 

(CAL/CMXX2-S-K) 


0 . 1000E + 01 
0 . 1000E+01 


0 .OOOOE+OO 
0 .OOOOE + OO 


1 . oooo 

1 .0000 


0 . 240E + 00 
0 . 2AOE+00 
0 . 2A0E + 0Q 


1 .0000 
1 .0000 


O.10QE-O3 0.600E-04 
0.ia0E-03 0.600E-04 
0.10UE-03 0.600E-0A 

HTACOR TAW 
(K) 


1 .0000 
1 . 0000 


0 . 1 602E+09 
0 . 3063E+03 


EXAMPLE PROBLEM 2 HEAT GENERATION - SI UNITS 


SI UNITS 

USER SPECIFIES PHYSICAL PROPERTIES (G, CP, MU AND K) 

COOLANT FLOW RATE IS SPECIFIED 

ROTATIONAL SPEED IS O.OOOE+OO RPM 

OUTPUT IN ORIGINAL NODE /INTERVAL FRAME OF REFERENCE 


Y 

(CM) 

0 .OOOOE+OO 
0 . OOOOE+OO 
0 . OOOOE + OO 


INPUT SLICES 

SLICE SPECIFIED HEAT 

NODE IN 

NO. IN 

INTERVAL 

NO 

. INTERVAL 


( J/KG) 

1 

125 

125 0 . 

13293E+07 

2 

125 

250 -0. 

1 3293E + 07 

I 

OPTIONS (IOPT) Z 

R 


FI HT FF 

CR (CM) (CM) 



0. OOOOE+OO O.OOOOE+l 

1 

0 0 0 

0 0 . 1GQQE + Q2 O.OOOOE + l 

2 

0 0 0 

0 0.2000E+02 O.OOOOE+l 

I 

T 

TT 

P 


(K) 

(K) 

(M/CM**2) 


0 . 2976E+03 

0 . 3000E+03 

0 . 9725E + 02 1 

1 

0 . 1A61E + 0A 

0 . 1 625E + 09 

0 . 57476+02 1 

2 

0 . 2975E+03 

0 . 2999E + 03 

0.9724E+02 l 


RE 

PR 

ST 

1 

0 . 1561E+07 

0 .719AE + 00 

0. OOOOE+OO l 

2 

0 . 1561E + 07 

0.719AE+00 

0. OOOOE+OO 1 


A 

(CMXK2) 

0 . 1 000E + 02 
0 . 1 000E + 02 
0. 1000E + 02 


DH 

(CM) 

0 . 3568E+ 0 I 
0.3560E+01 
0 . 3568E + 0 1 


TW 

(K) 


0 . 5330E+ 03 
0 . 533OE+03 


DX 

(CM) 


0. 1000E + 02 
0 . 9999E + 01 


XTOT 

(CM) 


o . inonE»02 
0 , 20QUE + 02 


0 . OOOOE+OO 
0. OOOOE + OO 125 
0 . 0000E * 00 125 


PT 

(N/CMX*2) 

0 . 1000E+03 
0 .8399E+02 
0 .9999E+02 


0 . OOOOE + OO 
0 .OOOOE+OO 


W 

(KG/S) 

0 .787AE+00 
0 .7879E + 00 
0 . 7874E+O0 


0 . OOOOE+OO 
0 .OOOOE+OO 


V 

(M/S) 

0 .6915E + 02 
0 .57AAE + 03 
0 .6914E + 02 


GAMMA CP 

CALZ(GR-K) 


vise . 

GR/ ( CM-S ) 


TH. CNDT. 
CAL/(CM-S-K) 


0 . 2Q00E + 00 
0.7497E+00 
0 , 2000E+00 


0. IA0E + 01 
0. 140E + 01 
0. 1A0E+01 


FCOR 


NC 

(CAL/CMX*2-S-K) 


0. 1 000E+01 
0. 1 OOOE + Q1 


0 . OOOOE + OO 
0 .OOOOE + OO 


1 .0000 

1 .0000 


0.24QE+00 
0 .2A0E + 00 
0 . 2A0E+00 


1 . OOOO 
1 . OOOO 


0.1806-03 O.6OGE-04 
0.100E-03 Q.600E-0A 
0.160E-03 0.600E-G4 

HTACOR TAW 
(K) 


1 . OOOO 

1 . oooo 


0 . 1603L+OA 
0 . 3049E + 03 
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